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CHAPTER  I 

INTRODUCTION 

We  have  attempted  to  demonstrate  the  Wigner  cusp  phenomenon 
in  neutron  elastic  scattering  and  total  cross  sections  near  the  thresh- 
olds for  excitation  of  the  first  excited  states  of  the  dominant  iso- 
topes of  Li,  Fe,  Zr,  Ba,  and  Ce .   The  energy  region  investigated 
extended  approximately  300  kev  on  either  side  of  threshold. 

Theoretical  work  (see,  for  example,  Newton,   Baz,  '   Fonda, 

6,7x 
and  Meyerhof   )   provides  detailed  predictions  of  possible  cusp  effects. 

The  effects  originate  essentially  in  two  fundamental  physical  principles; 
the  conservation  of  incident  flux,  expressed  in  the  theory  as  unitarity 
of  the  scattering  matrix ;  and  the  uncertainty  principle  of  Heisenberg, 
expressed  in  the  theory  as  analyticity  of  the  scattering  matrix.   The 
magnitude  and  shape  of  cusp  anomalies  can  yield  information  useful  in 
scattering  analysis  and  in  establishing  spins  and  parities  of  excited 
states . 

Cusp  anomalies  have  been  seen  in  proton  elastic  scattering  near 
thresholds  for  (p,n)  reactions.  '  '    The  theory  shows  that  one  essen- 
tial requirement  for  a  sharp  cusp  is  that  coulomb  effects  be  absent 
from  that  reaction  channel  which  provides  the  threshold.    Any  reaction 
with  outgoing  neutrons,  such  as  (p,n)  (d,n)  (a,n)  (n,nT),  satisfies 
this  requirement 


The  advantage  of  studying  Wigner  cusps  in  elastic  neutron, 
rather  than  elastic  proton  scattering,  for  example,  is  the  easier 
analysis  caused  by  the  absence  of  coulomb  effects  in  the  incoming 
channel.  After  estimating  possible  cusp  magnitudes,  we  concluded  that 
the  effects  of  inelastic  neutron  scattering  reactions  on  elastic  neu- 
tron scattering  might  be  quite  noticeable  for  selected  cases.  We 
hoped  thereby  to  find  a  method  for  spin  assignment  of  excited  states 
and  to  elucidate  interaction  mechanisms  between  nuclei  and  neutrons 
of  few  Mev  energy. 


CHAPTER  II 

THEORY  OF  WIGNER  CUSPS 

Previous  articles  on  the  theory  of  Wigner  cusp  phenomena  have 
generally  emphasized  the  simple  spin  cases  (e.g.,  spinless  particles); 
also,  the  results  are  not  always  given  in  a  form  readily  compared  with 
experiment.  We  would  like  to  supplement  the  previous  results  by  con- 
sidering more  general  spin  and  angular  momentum  configurations.   Our 
results  will  provide  a  convenient  form  for  experimental  comparisons. 
In  addition,  the  form  of  the  results  for  pure  Wigner  cusp  effects  will 
be  particularly  suitable  for  subsequent  deduction  of  "energy-averaged 
cusp  expressions,  as  first  derived  in  Ref.  6. 

Pure  Wigner  Cusp  Theory 

We  will  compute  the  effect  of  reactions  of  the  type  X(n,n')X* 

on  neutron  differential  elastic  and  total  cross  sections  near  threshold 

12 
energies.   Blatt  and  Biedenharn   have  derived  necessary  general  ex- 
pressions for  cross  sections.   For  completeness,  we  give  in  Appendix  A 
a  summary  of  the  physical  basis  for  the  cross  section  formulas,  start- 
ing with  the  wave  functions  and  the  scattering  matrix,  S . 

In  applying  these  formulas  we  assume  that  elastic  scattering 
with  change  of  orbital  angular  momentum  Z     may  be  neglected.   Because 
of  the  parity  rule  requiring  M   =   0,  2,  h    .  .  .,   elastic  scattering 
with  change  of  £     turns  out  to  be  an  unimportant  phenomenon  in  most 


13 
situations.    In  particular,  in  the  present  experiment  this  restric- 
tion is  very  mild,  since,  with  one  exception,  the  target  nuclei  were 
even-even.  Hence,  with  I  =  0,  i  =  \,    s  =  \   (i,  i,  and  s  are  target, 
projectile,  and  channex  spin,  respectively),  there  is  no  possibility 
of  elastic  scattering  with  change  of  £.      The  exceptional  case  is  that 
of  Li   (I  =  3/2).   In  this  case,  NL   ^  2  is  limited  by  penetrability 
effects.   For  example,  although  incident  £  =  0     waves  theoretically 
can  be  scattered  elastically  into  £   =   2  outgoing  waves  (with  com- 
pound spin  J  =l),  it  happens  that  in  the  energy  region  of  interest 
kR  =  0.3  an(i  any  process  involving  £  =   2  is  highly  improbable.  Con- 
sequently, for  the  present  experiment,  specializing  the  derivation  to 
A£  =   0  is  not  a  restriction  on  the  validity  of  the  results. 

If  one  considers  angular  momenta  situations  in  which  there  can 
be  scattering  involving  change  of  channel  spin  s,  but  no  change  in 
channel  energy,  complexities  arise  in  a  derivation  of  cusp  effects. 
To  preserve  clarity  in  the  presentation,  we  will  restrict  our  deriva- 
tion to  situations  in  which  channel  spin  s  is  separately  conserved, 
in  addition  to  the  previous  assumption  that  orbital  angular  momentum 
i  is  conserved.   It  will  be  shown  that  this  restriction  does  not  af- 
fect the  validity  of  the  derivation  as  it  applies  to  present  experi- 
mental results.  We  will  indicate  the  nature  of  the  additional  complex- 
ities in  cusp  theory,  which  are  present  without  this  restriction,  at 
a  suitable  point  in  the  derivation. 

The  following  equations  are  numbered  3.l6,  k. 5,  and  k.6   in 
Ref.  12.   We  have  specialized  them  to  elastic  scattering  with  no  change 


of  a,    l,    s  or  i .  Accordingly,  the  differential  elastic  scattering 
cross  section  is  given  by 

I+i 
da  „  =  Z     4  (2s  +  1)/[(2I  +  l)(2i  +  l)]|-  da (l) 


j         =  E     j(2s  +  1)/[(2I  +  l)(2i  +  1)][  da 
a'a      b=  I-i  L  J     ' 


as 


■where 


?  °° 

da  =   |>     /(2s  +  1)]      Z     BT(as;as)   PT(cos  Q)dSl  (2) 

OCS ;  OtS  (X  j  _(-.      ij  -L" 


B  (QBjas)    =  (l/U)      S        £        £       E 

Jl     J2      'l        *2 

X  Z2(i1  Jx  i2  Jg,    s   L)  (3) 

x  Ret(1  "  SaLl5c^/  (1  "  8«igjaBi2» 

where  a  refers  to  the  type  of  incident  particle  and  the  state  of  the 
struck  nucleus;  -ft   is  (2it)   X  the  deBroglie  reduced  -wavelength  of 
the  incident  neutron/target  nucleon  system;  s,  I,  and  i  are  channel, 
target,  and  projectile  spins,  respectively;  and  I,    J  are  orbital  and 
total  angular  momentum,  respectively.  Subscripts  1  and  2  are  dummy 
indices  for  use  in  summation.  All  sums  are  unrestricted,  except  that, 
in  practice,  only  one  such  sum,  for  example,  J, ,  runs  to  infinity,  be- 
cause of  the  vanishing  of  Z  coefficients.   Re  stands  for  real  part 

of  the  bracketed  expression.   S   .    „  is  the  S-matrix  element  for 

*  asi;asi 

elastic  scattering  with  £     and  s  conserved  ( see  Appendix  A) . 


L-V-^p         1111 
Z(ix  Jx  ^  J2,  s  L)  =  i   X  d(2ix  +  l)2  (2i2  +  l)2  (2JX  +  1)2(2J2  +  l)2 


X  W(ix  Jx  Z2   J2,  s  L)(i1  i2  0  Ol^  ig  L  0) 

lU  15 

where  W  is  the  Racah  coefficient.    Tables  of  Z  are  available » 

16 
The  properties  of  W  are  discussed  by  Biedenharn,  Blatt,  and  Rose. 

The  energy  dependence  of  the  elastic  cross  section  is  contained 
in  S.  We  will  use  the  unitary  and  analytic  properties  of  the  S-matrix 
to  obtain  an  expression  for  the  elastic  cross  section  in  the  vicinity 
of  a  reaction  threshold. 

Among  the  inelastic  reactions  which  can  proceed  from  a  given 
incident  partial  wave,  we  consider  only  the  one  with  lowest  threshold 
energy.  Note,  however,  that  an  excitation  of  a  higher  excited  state 
is  not  excluded  as  long  as  its  threshold  has  the  lowest  energy  of  those 

reactions  which  may  proceed  from  a  particular  incident  partial  wave. 

3 

Baz  showed  that  the  procedure  for  including  higher  inelastic  thresh- 
olds of  the  same  dominant  partial  wave  is  straight-forward  in  principle, 
although  much  more  complicated.   Since  the  present  experiments  were  done 
in  energy  regions  near  first  excited  states,  the  case  we  consider  is 
adequate,  and  the  physical  content  of  the  results  is  not  obscured  by 
these  complications . 

Also,  other  true  inelastic  processes  (e.g.,  neutron  capture) 
occurring  at  lower  energy  than  the  threshold  energy  of  interest  are 
neglected. 

In  addition,  we  restrict  the  energy  region  in  which  the  deriva- 
tion is  valid  to  that  in  which  the  reaction  is  dominated  by  a  single 


incident  partial  wave.   In  practice,  such  an  energy  region  is  usually 
sufficiently  large  to  be  of  experimental  interest.  For  example,  an 
incident  d  wave  contributes  more  than  90  percent  of  the  cross  section 
over  an  energy  interval  of  150  kev  above  threshold  for  inelastic  neu- 

+  1U0  17 

tron  scattering  leading  to  the  2  first  excited  state  of  Ce    .    One 
can  explain  this  in  terms  of  the  optical  model  by  noting  that  Ce  has 
a  large  transmission  coefficient  for  (low  energy)  outgoing  s   wave 

-1  o 

neutrons,   and  that  sufficiently  close  to  threshold,  inelastic  neutron 
scattering  is  always  dominated  by  that  incident  partial  wave  which  is 
associated  (through  angular  momentum  conservation)  with  an  outgoing 
s   wave . 

Finally,  to  further  simplify  the  presentation,  but  not  as  an 
essential  limitation  on  the  theory,  we  will  only  consider  reactions 
dominated  by  outgoing  _s  waves.   The  complete  theory  in  this  respect 
shows  that  reactions  dominated  near  threshold  by  outgoing  waves  of  or- 
bital angular  momentum  Z   ,   greater  than  zero,  produce  cusp  effects 

in  the  (i  )th  derivative  of  the  elastic  cross  section.    Therefore, 

o'  ' 

these  cases  have  less  experimental  interest.  With  one  exception,  this 
restriction  is  not  essential  since  we  measured  threshold  effects  of 
reactions  with  outgoing  s   waves.   The  exceptional  case  is  the  reaction 

Zr  (n,n')  Zr    (1.75  Mev)  which  Tucker,  et  al.  have  found  dominated 

19 

by  outgoing  p  waves  already  at  about  30  kev  above  threshold.    The 

90 
necessary  modifications  to  the  theory  in  the  case  of  Zr   will  be  ob- 
vious . 
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Under  the  above  assumptions  the  reaction  [see  Eq.  (5)  below] 
may  be  described  by  its  partial  reaction  cross  section 

where  primes  identify  the  particular  channel  quantities  which  couple 

to  i',   the  dominant  incident  partial  wave  and  g  ,  =  (2J1  +  l)/(2s'  +  l) 

s 

is  a  statistical  factor  obtained  after  coupling  angular  momenta,  averag- 
ing over  initial  magnetic  quantum  numbers  and  summing  over  final  magnetic 
quantum  numbers  (see  Blatt  and  Weisskopf   or  Blatt  and  Biedenharn   ) „ 

To  simplify  the  notation,  we  now  drop  a  number  of  subscripts. 
First  a'   is  superfluous  since  we  always  consider  a  neutron  incident 
upon  the  ground  state  of  a  target  nucleus.   Second,  we  reduce  the  double 
subscript  set  on  the  S-matrix  element  to  a  single  set;  the  meaning  of 
S,   as  given  in  Appendix  A,  should  be  kept  in  mind. 

Equation  (4a)  becomes" 

Equation  (4b)  is  essentially  the  statement  that  the  scattering  matrix 
J3  is  uni'cary  (see  Appendix  A).   The  corresponding  expression  for  the 
total  reaction  cross  section  is; 

1+1  2s'  +  1      i,+S'       J' 

°r   "  *ri«  -  §  =M|  (21  ?  l)(2i  +  l)  ^.^  <JV   <*> 

Inverting  Eq„  (4b),  the  S-matrix  element  may  be  written; 

#,,.(«* 21 5j:,,xi- ,£...>*  (6) 


Jl 

Equation  (6)  defines  26.  t  t  as  the  (real)  phase  of  the  S- 

Jj    s 

J' 
matrix  element  S.,  ,  .   To  simplify  the  notation,  we  have  expressed 
jL    s 

J*  2  J1 

a   i  I   in  units  of  it  k     a    ,  .   This  normalization  will  be  used  in 
r£   s  &s 

the  remaining  discussion  except  in  instances  where  clarity  requires 
that  the  units  be  specifically  restored. 

Equation  (6)  is  certainly  valid  above  threshold.   The  usual 

2  3 

procedure    is  to  expand  the  complete  phase  (which  may  have  real  and 

j' 
imaginary  parts)  of  S  ,  t   in  energy  about  threshold.   One  consequence 

4/   S 

of  this  procedure  is  that  the  fixed  value  of  the  real  phase  at  thresh- 
old appears  in  the  final  expressions.   On.  the  other  hand,  Eq.  (6)  is 
quite  general  and  there  is  no  reason  to  assume  that  the  form  of  Eq.  (6) 
should  not  persist  even  away  from  threshold  (above  threshold).   Conse- 
quently, we  shall  allow  the  real  part  of  the  phase  of  S,   i.e.,   26, 
as  well  as  a  ,      to  vary  over  the  energy  region  of  interest,  if  the 
physical  situation  requires  that  they  do  so.   This  leaves  us  free  later 
to  examine  the  consequences  of  energy  averaging  over  a  cusp  region  where 
the  phase  26  and  the  cross  section  0  might  be  fluctuating  rapidly 

because  of  the  presence  of  many  narrow  resonances. 

7 
Using  R-matrix  theory,  Meyerhof  has  shown  that  the  S-matrix 

can  always  be  written  in  a  form  equivalent  to  Eq.  (6)  if  £     =0,  inde- 
pendently whether  the  energy  is  above  or  below  threshold.  An  analytic 
continuation  in  energy  across  threshold  always  exists,  even  in  the 
presence  of  many  narrow  and  possibly  overlapping  resonances  which  would 
produce  rapid  variation  with  energy  in  the  phase  and  in  the  reaction 
cross  section. 
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The  prescription  for  continuing  Eq.  (6)  below  threshold,  assum- 
ing outgoing  s   wave  neutrons  in  the  reaction,  is  given  by  Eq.  (29a)  in 
Ref„  7;   essentially,  one  may  let  a     =   \o    \      above  threshold  go  to 
i | a  |   below  threshold.   It  is  interesting  that  in  R-matrix  theory  this 
prescription  actually  results  from  a  change  In  the  effective  logarithmic 
derivative  of  the  reaction  channel  wave  function  from  a  pure  imaginary 
number  above  threshold  to  a  pure  real  number  below  threshold.  This  is 
related  to  the  fact  that  outgoing  wave  function  (exp  ik  r)/r  has  to 
change  to  (exp  ~|k  |r)/r  as  the  energy  (E)  is  changed  from  above 

threshold  (Em)  to  below.   k  =  n/2M(E  -  Em)/ft  ,   M  =  reduced  mass 
T  o  T  '        * 

of  tne  outgoing  neutron.   This  prescription  becomes  more  complicated  if 

one  considers  outgoing  waves  in  the  reaction  channel  with  £     £   0.   The 

o  ' 

necessary  modifications  for  extension  of  the  present  theory  to  these 

cases  is  also  given  in  Ref .  7* 

Using  the  fact  that  a  . ,  ,  <  <  1  near  threshold,  we  can 

ri  s  •  ' 

simplify  Eq.  (6)  to 

and,  upon  substituting  the  analyticity  prescription  into  Eq.  (7),  we 
obtain 


BSv-<«**»iv>u-i<£v0}> 


(8) 


where,  following  Newton,  -j    Y    means  A  if  E  >  E  ,   and  B 

E  <  E  ,   where  E    is  the  reaction  threshold  energy.   It  should  be 
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noted  that  a      ,  ,   in  Eq.  (8)  is  strictly  equal  to   a   ,  f  (  | E  -  E   |) 

(  l\  ^ 

the  accompanying  ■<    .  r     contains  the  analytic  properties  of  the  expres- 
sion. 

j» 
We  have  separated  S.,  ,   into  two  terms .   The  first, 

Ju     s 

exp  2i  8,1    i    ,      can  he  interpreted  as  the    "S-matrix  element  -without 

iu    S 

•7 

reaction  effects,"  if  i  =0.   When  this  term  is  substituted  into 

o 

Eqs .  (l),  (2),  and  (3)>  and  the  indicated  sums  are  performed,  the  re- 
sult is  the  differential  elastic  cross  section  without  reaction  ef- 
fects, "  which  we  designate  by  da    .  The  second  term, 

j«  ji         f  l  I 

-(exp  2i  6»i    <)   2   a    of    »    l    •     f      produces  a  cusp  term  after  substitution 
x>   s  r  xj   s      ill 

into  the   same  equations. 

Accordingly,    the  Wigner  cusp  term     A  da^  >/v     is  defined  by  the 
set  of  equations 


a;  a 


da         =  da°       +  A  da  (Q) 


where 


I+i 
AdVa  =       Z,         ,[(2s'   +  1)/(2I  +  l)(2i   +  1)]  A  da      ,         ,  (10) 

S     =    1-1  ' 


2 

Adaas'-as'  =  C*a  ^2s'  +  l)]     z    AB^as'jas')  PL(cos  e)dfl      (11) 

J_i=U 
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ABT(as';as')  =  (l/k)  \    Z  Z  Z  Z2(in  J.  iv  J',  s'  L) 

L  l.T   «  ,T' 


J1J1J 


X  Re 


1  T '      IT1 

(1  -  exp  21  6i  g,)*  (exp  2i  5,,^,)  5-  cr^^, 


i's"  2  ri's'   i 


+  Z  Z  Z  Z2(i'  J'  £     J  ,  s'  L) 
J'  J2  ^ 


X  Re 


^(exp  2i  6iIsl)  ^  «£jV  |ij  (1  ~  ^  2i  8    ,)J  j 


J1  2  J1 

where  a  „t  1   is  still  expressed  in  units  of  it  -ft  g  ,   [see  Eq.  (6)] 
r  *  s  s 

The  sums  over  £   ,    J,,  £   ,   J   in  Eq.  (3)  are  reduced  "because  of  the 
assumption  that  a  single  wave  £%      is  dominant  in  the  reaction. 

Because  Z(£     J  i  J-,  s  L)  =  (-l)L  z(  ig  Jg  £,  J,,  s  L)  and 
Re(A"x"  B)  =  Re(A  B*),  we  may  combine  these  sums,  changing  dummy  indices 
and  dropping  the  extraneous  index  (X     for  simplicity,  to  obtain 


A  BT(s';s')    =  (1/10    Z  Z  Z     2   Z   (i  J   i'    J*,    s1    L) 

L  JIJ1 


X  Re 


(1  -   exp  2*  6^,)*   (exp  21  5^,    ,)  g  D^.    .  -j* 


(l/U)   Z  Z  Z     Z  (i  J  i'    J',    s!    L) 

JIJ1 


J  J' 

X  2   sin  5  „    ,    a    . ,    ,   Re 
is        r£  s 


<£}■«»  1   (SbJI    ,   -8J    ,) 


Taking  the  Re  [  ]   in  the  last  equation  yields 
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A  BT  ( s  ' ,  s  ' )  =  -  ( l/k)   EEZ  Z2(  £   J  £ '  J ' ,  s  '  L) 
L  J  i  J* 


jsin(26^s,  -  6^, 
X  a  .,  ,2  sin  6fl  ,  i      Ti      T   f  ^J 


-os(2Vs'  "  5is' 


Within  the  limitations  imposed  during  the  derivation,  Eqs.  (9), 
(10),  (ll),  and  (12)  describe  expected  Wigner  cusp  effects.   The  cusp 
effect  is  completely  contained  in  the  coefficients  BT   in  Eq.  (2); 
that  is,  B  =  B°  +  A  B  ,   where  A  B   is  given  by  Eq.  (12). 

Since  it  is  frequent  practice  (cf .  Ref .  20 )  to  analyze  neutron 
scattering  data  in  terms  of  the  BT   coefficients,  the  utility  of  Eq. 
(12)  for  cusp  interpretation  is  clear.  Should  a  cusp  effect  be  found 
in  only  one  or  two  of  several  experimental  B   coefficients  then  the 
properties  of  the  Z  coefficients  severely  restrict  the  participating 
partial  waves . 

We  now  recall  that  the  cusp  predictions,  Eqs.  (9)  to  (12)  are 
not  generally  valid  if  scattering  is  possible  which  conserves  energy 
but  not  channel  spin,  because  most  detection  systems  are  only  energy- 
selective.   In  such  cases,  the  threshold  reaction  is  coupled  (via 
unitarity  of  the  S-matrix)  to  other  off-diagonal  (different  initial  and 
final  state  channel  spin)  elements,  as  well  as  a  diagonal  element 
("pure"  elastic  scattering)  of  the  S-matrix  (see  Eq.  (ha)).     An  addi- 
tional sum  over  final  state  channel  spins  is  required  in  Eqs.  (l)  to 
(3).  As  a  result  additional  cusp  terms  appear  in  Eq.  (ll).  Unless 
polarization  is  also  measured,  comparisons  of  theory  and  experiment 


Ik 


are  ambiguous  in  such  cases .  None  of  the  present  experimental  results 
required  this  extension  to  cusp  theory  for  comparison  purposes. 

The  expected  Wigner  cusp  in  total  cross  sections  (a  )  is 
easily  obtained  if  we  recall  that,  by  assumption,  the  interaction  in 
the  (i')    incident  channel  is  purely  elastic  below  threshold. 


aT  =  /  (da/df2)dfi  +  ori,  |q| 


/  (A  da/dH)dfi  +  ari,  |q 


where  a  °  is  the  total  cross  section  "without  reaction  effects." 

Using  the  fact  that     /   P  (cos  0)dft  =  6?"  _   kit      (where     5?  _ 

■Li  Li  j  0  L,  0 

is  a  Kronecker  delta)  we  obtain 

*"*  ~s-fii.il  (si '2i)(2iH  i)  (2r^Ty^ABo(s,^')+°rr{o}  ^ 

Substituting  Z(i  J  V   J',  s'  0)  =  6^,  6jjt(-l)J"S'  (2J  +  l)?,  where 

K  K 

5--i  and  6   ,   are  Kronecker  deltas,  into  the  equation  for  AB 

(with  L  =  0),   Eq.  (12),  and  putting  this  result  into  Eq.  (ik) ,   one 

obtains 


J?        (2b'+1)      ^'+S'    At*2(2j'fl)   J' 
T  =  ..t|T  -I  (21  +  l)(2i  +  1)   Z,      A  (2s,+l)    ari's' 


■'"^I'vli 

-  sin25rs, 
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where  we  have  explicitly  indicated  the  units  of  a    .,    ,   (see  Eq.  (6)) 

I*  Jj    s 

Finally,  recognizing  that  the  summations  in  Eqs.  (15)  and  (5)  are 
identical,  we  may  combine 


A 

r2  sin26k'j 

r  cos25iV 

oJ  ' 

^sin26lV 

I  -   sin  25^ 

s' 

and  obtain  the  result 

^        (2s'  +  1)       i'+S'    f*  fc2(2J'  +  1)'  J 
T     .'!|I-1|  (21  +  l)(2i  *   1)  j«!|i..s'|V  (Sb'+I)   °r*V 


cos  2&J,s, 


-  sin  26  „ , 


.J' 
Vs 


(16) 


Equation  (l6),  combined  with  Eq.  (13)>  is  the  general  expression  for 
Wigner  cusp  effect  in  the  total  cross  section,  within  the  limitations 
imposed  by  the  derivation. 

Wigner  Cusps--Spinless  Particles 

We  now  extract  the  special  case  of  spinless  particles  from  these 
expressions,  for  comparison  with  other  work. 

If  all  spins  are  zero,  including  that  of  the  residual  nuclei 
in  the  reaction,  then  Eq»  (12)  becomes 


2  0   fs 

A  BL(0;0)  =  -.(1/4)  2  Z  (4  16   0  0,  0  L)  2  sin  b£  a^  -I 


sin(250  -  8,, 


os(260  -  6 
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LP,  - 

Since  Z(i  I  0  0,  0  L)  =  (-1)    (21  +  l)2  (l  0  0  0   i  0  L  0) 

K  K 

=  6  .,T  ,  vhere  5    is  a  Kronecker  delta, 

the  sum  in  A  BT   is  reduced  to  one  term,  and  using  Eqs .  (10)  and  (ll) 

Li 

we  obtain 

2  ,         "   0  f  sin(26o  "  5P  1 

Ada=-(*  A)  S  a^  2  Bins,  |  }  PJ  drt 

i=0  ^  cos(26Q  -  6.)  J 

where  6.  =  6.-.   This  expression  can  be  written  as 


A  da  s  2  Re 


thus 


/  0         fc2x 

(i  1  *    ^  (1  -  exp  tt  B/)p/       g;  |  exp  21  5q  {^ 


A  da  =  2  Re 


(lA)**(9)i«W2i80{j}(ar) 


(17) 


where  f*(0)  is  the  complex  conjugate  of  the  scattered  amplitude,  with 

its  explicit  form  evident  by  comparison  with  the  preceding  line,  and 

2   O 
where  (a)    =  («  k  a  _);  i.e.,  the  experimental  reaction  cross 
r  exp         rO  ' 


section, 


k 

This  is  just  the  cusp  term  obtained  by  Baz,   for  the  case  of 


spinless  particles.   Baz  emphasized  the  explicit  appearance  of  the  am- 
plitude, suggesting  that  there  is  additional  theoretical  significance 
in  cusp  data  beyond  cross  section  measurements  (squared  amplitudes) 
that  are  usually  made .  We  would  like  to  caution  that  most  cases  in 
practice  involve  non-zero  spins;  consequently  scattered  amplitudes, 
in  general,  are  azimuth  angle  (cp)  dependent.  Comparison  of  theoretical 
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cusp  effects  with  experiment  may  be  undertaken  if  the  theoretical  ex- 
pressions have  been  averaged  and  summed  over  magnetic  quantum  numbers. 
This  is  essentially  what  we  have  accomplished,  and,  in  general,  little 
vesitge  of  complete  amplitude  forms  is  found  in  the  resulting  expres- 
sions . 

On  the  other  hand,  it  is  useful  to  see  from  Eq.  (17)  that  one 
can  estimate  expected  Wigner  cusp  magnitudes  from  the  expression 

i_ 
|  A  dcx|  =  a  (da)2/\ 

the  remainder  of  the  cusp  expression  being  of  order  unity.   The  ex- 

l/2 
pected  percentage  cusp  effect  is  inversely  proportional  to  (da)  ' 

and  exploratory  experiments  might  well  be  undertaken  at  angles  where 

the  cross  section  is  least. 

The  spinless  particle  result  for  total  cross  sections  is  a 

trivial  reduction  of  Eq.  (l6)  to 

f  cos  2&n  1 

T    T      r'exp  \    .   OR  J 

**•  -  sin  28nJ 

2  6 
This  result  is  in  agreement  with  those  of  others.  ' 

This  concludes  our  discussion  of  pure  Wigner  cusp  phenomena. 
Two  other  special  cases,  applicable  to  the  present  experiment,  are  de- 
duced in  Appendix  Ik 
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Theory  of  "Energy-Averaged"  Wigner  Cusps 

We  turn  now  to  the  problem  of  calculating  the  expected  threshold 
effect,  if  the  measured  cross  sections  are  effectively  energy  averages 
of  the  actual  detailed  cross  sections,  containing  closely  spaced  reson- 
ance structure.   That  the  resulting  effects  are  related  intimately  to 

6,7 

pure  Wigner  cusp  phenomena  was  first  discussed  by  Meyerhof .,     He 

showed  that  the  energy  average  of  Wigner  cusp  expressions  leads  to  a 
predicted  threshold  effect  which  we  shall  call  an  "energy-averaged" 
Wigner  cusp. 

The  physical  situation  we  are  considering  is  that  the  reaction 
and  scattering  cross  sections  are  dominated  by  many  closely  spaced 
resonances .   By  closely  spaced  we  mean  a  large  number  within  the  energy 
interval  A  E  over  which  the  averaging  is  performed:   either  experi- 
mentally, because  of  finite  incident  beam  resolution  or  because  of 

overlap  of  the  resonances;  or  theoretically,  by  choosing  A  E  >  >  D, 

7 
where  D  is  the  average  level  spacing.  No  essential  restrictions 

will  be  placed  on  the  level  width  to  spacing  ratio,   r/D. 

Noting  that  the  energy  dependence  of  differential  elastic  cusp 

effects  is  contained  in  ABJs'js1)  (see  Eq.  (12))  we  can  isolate  the 

energy  dependent  quantities: 

T    r    Sln(25i's'  "  ^l 

<tfr*V2si*5is<{         rpRJ»  RJ     A)  ™ 

Lcos(26i?s,  -  ois.) 

Brackets,   (  ),   imply  that  an  energy  average  is  to  be  taken 
over  an  interval  containing  many  resonances.   The  phases  8    and 

i/S 
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ji 

5 . :  i   are  rapidly  varying  with  energy  over  the  averaging  interval  and 
Xf    s 

we  will  assume  they  are  uncorrelated  for  J  f   J1   and/or  Z  f  &' .     This 
is  the  well  known  statistical  assumption  for  the  compound  nucleus.   In- 
tuition might  suggest  that  uncorrelated  phases  would  lead  immediately 
to  breakup  of  expression  (l8)  into  a  simple  form,  for  example,  the 

product  of  the  averages  of  the  factors ,   It  turns  out  that  such  a  break - 

j« 

up  is  by  no  means  obvious,  because  of  the  presence  of  a    . ,    ,  within 

r  Jo   s 

the  averaging  brackets . 

Using  R-matrix  theory  and  applying  the  random  sign  assumption 

7 

to  certain  value  quantities  y  ,   which  will  be  discussed,  Meyerhof 

A, 


showed  that 


2     j '  j  • 

-2  sin     6,i_  n  f   (    oj„ , 


TJi    s 


,ln2£Vs' 


( £..,        T?v  >  -       rrs    y  (19) 


and 


sln  26i's' 
The  value  quantities  y   mentioned  above  were  introduced  by 

A, 

21 
Wigner  and  Eisenbud   to  parametrize  the  terms  of  the  R-matrix,  each 

of  which  has  a  resonance  denominator.   Thus,  these  quantities  corre- 
spond to  resonances  and  could  be  indexed  in  our  notation  by  T  .  , 
where  X.  specifies  a  particular  resonance  of  the  set  associated  with 
angular  momenta  J,  i  and  spin  s.   The  random  sign  assumption  leading 
to  Eqs.  (19)  and  (20)  is  made  for  value  quantities  of  different  index 
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J,  i,   or  s,   but.  is  not  justified  for  relating  value  quantities  of 
the  same  J,  £     and  s,   but  different  X;   i.e.,  different  resonances 
initiated  by  the  same  partial  wave  may  have  correlated  phases .   In 
particular,  the  reaction  cross  section  appearing  in  the  averaging 
brackets  in  Eqs .  (19)  and  (20)  proceeds  through  resonances  correlated 
in  phase  to  the  remaining  expression  in  the  bracket,  and  it  is  surpris- 
ing  that  the  simple  results,  Eqs.  (19)  and  (20)  were  obtained. 

In  Appendix  C,  we  show  how  Eqs.  (19)  and  (20)  can  be  applied 
to  expression  (l8)  to  give  the  result 

/J'  0      .      RJ        fsin(2&^',    ,    -  bJ0    ,)\  v 

/      J'  \ 

ri   S  bKaol    bKTTl  (21) 


££'      JJ 

where  5   is  a  Kronecker  delta. 

Part  of  the  proof  of  Eqs*  (19)  and  (20),  which  led  to  Eq.  (21), 
depends  upon  a  cancellation  of  path  integrals  used  in  the  energy  averag- 
ing process,  which  will  be  complete  only  if  normally  slowly  varying 
quantities  remain  slowly  varying  throughout  the  averaging  interval. 
For  the  case  of  outgoing  s   waves  in  the  reaction  channel,  one  of  these 
quantities  (the  logarithmic  derivative  of  R-matrix  theory)  does  not 
vary  slowly  right  at  reaction  threshold.   In  essence,  this  difficulty 
is  just  the  result  of  rediscovering  pure  Wigner  cusp  phenomena  as  one 
formally  approaches  very  close  to  threshold  in  the  averaging  process. 
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This  means  that  the  derivation  of  Eqs .  (19)  and  (20)  is  not  strictly- 
valid  arbitrarily  close  to  threshold.  The  actual  experimental  result 
would  depend  on  the  exact  energy  distribution  of  the  beam.  As  long  as 
the  energy  distribution  of  the  beam  does  not  overlap  the  threshold, 
the  results,  i.e.,  Eqs.  (19)  and  (20)  and  consequently  Eq.  (21),  will 
be  applicable . 

After  justifying  Eq.  (21)  we  have  completed  the  derivation  of 
"energy -aver aged"  cusp  expressions.   Recalling  Eq.  (12),  which  contains 
expression  (l8),  and  using  Eq.  (21),  we  obtain 

(aJ\   ) 
(  A  BT(s';s')  )  =  -  (1/4)  Z     Z2(i'  J'  r    J',  s'  L)  {     riV    \   (22) 
L  J'  L   0     -I 

K  K 

the  sums  over  i  and  J  have  vanished  because  of  5-.,   and  6TTt 

in  Eq.  (2l).   The  energy-averaged"  differential  cross  section  cusp  term 
is  obtained  by  substituting  Eq<  (22)  in  Eqs.  (9-ll). 

Substitution  of  relation  (20)  in  Eq.  (l6)  gives  the  expected 
result  for  the  total  cross  section, 

0 


(AcT>  -Q  (23) 


That  is,  there  is  no  energy- averaged  cusp  in  the  total  cross  section. 

22 
This  is  consistent  with  the  basic  optical  model  assumption,   used  in 

relating  the  model  to  experiment,  that  the  energy-averaged  total  cross 

section,  being  linearly  dependent  upon  the  scattering  matrix  elements, 

is  smooth  in  energy  even  across  reaction  thresholds.  We  will  come  back 

to  this  point  when  summarizing  the  theory. 
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Energy- Ave raged  Cusps--Spinless  Particles 

As  an  example,  for  comparisons  to  previous  work,  consider  the 
special  case  of  spinless  bombarding  and  target  particles 

I=i=s'  =0j   £'  =  J' 

Equation  (22)  becomes 

(  ABL(0;0)  )  =  -  (l/U)  Z2(i*  V    I1  £',   0  L)  j   ri '°  |        (22a^ 

Using  Z2(i'  i1  i1  i',  0  L)  =  (2i'  +  l)2  U'  /'  0  0  |  i1  i'  L  0)2, 
indicating  explicitly  the  units  of  (a  .tp.)*   and  by  substituting  Eq. 
(22a)  into  Eqs  .  (9-11),  we  obtain 

(Ada   )  =  -  (2J1  +  l)(jt  *2(2i'  +  l)  (ai*  >A«)2  U'  i'  0  0  |  I'  I'  I  0)2 

x  PT  &Q, 

Li 

,-(2i'+l)«ari,>e^An)(Pi,)2{J,  (24) 

where   (a  .,)     means  units  have  been  restored  to  (a    „,)      and  experi- 
ri  exp  ri 

mental  values  may  be  substituted  for  this  quantity.  Equation  (2^-)  is 
the  same  as  the  formula,  Eq.  (l6)  of  Ref .  6. 

"Energy-Averaged"  Cusps--Spin  \   Incident  Particles 

We  give  the  results  for  a  second  special  case  of  Eq.  (22)  of 
particular  interest  in  comparison  with  the  present  experimental  data, 
where  the  incident  particle  is  a  neutron,  the  target  nucleus  is  even- 
even,  and  incident  d  waves  dominate  the  reaction. 
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I   1 


I  =  0,  i  =  i,   j'  =  £'   ±  \,  i1  =  2,   and  s1  = 


£'+i 


(a  M£'   =2))  .--.(UrJ.^>  expA*)   £   Rj 


eXp     J'=i'-t 


S  [Z2(2J'  2J',  i  L)/(2J'  +  IJIPlJo}^       <25) 


X 

L 

where  RT,  ■  (a  x>   /(a   ,) 

J      ryef  exp'   ri  exp 

Although  not  immediately  obvious,  Eq.  (25)  is  identical  to  Eq. 

(161)  in  Ref.  6.   To  compare,  one  must  identify  R.,  i  =  (i'+l)/(i'+l 

i/  "•*  p 

+i'r),   R  ,   i  =  i'  r/(i'  +  1  +  i'  r),  where  r  is  defined  and  used 
*  "  2 

in  Ref .  6. 

Since  Eq.  (25)  is  important  to  present  experimental  comparisons, 

it  will  be  useful  to  give  the  terms  in  detail.   It  turns  out  that  the 

cusp  prediction,  Eq.  (25),  is  quite  insensitive  to  the  ratio  RT)  ; 

recall  that  R  ,  was  defined  as  the  fraction  of  the  reaction  probabil- 
J 

ity  which  originates  from  compound  resonances  of  a  particular  J',   for 
example  i1  +  J.  We  make  the  reasonable  choice  R  ,  =  i(2J'+l)/(2i'+l) 
which  is  the  same  as  the  fraction  of  the  total  number  of  incident 
states  of  the  ( £ ' )    partial  wave  which  have  total  angular  momentum 
J1.  With  this  choice  of  RTi,  Eq.  (25)  becomes 


(AMf=2))=-((,,Ti,JeJk«){$ 


X  (2J1  +  l)"1  Z  |  Z2(2J*  2J',  ^  L)  P_  dft 
J',L2  2     L 


X  [1  +  38/35  P2(cos  e)  +  18/35  P^(cos  e)]dfl      (26) 
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Equation  (26)  actually  applies  to  any  scattering  experiments 
on  medium  and  high  atomic  numbered  even-even  targets  at  energies  where 
random  phase  concepts  are  useful,  and  near  the  threshold  of  a  reaction 
dominated  by  incident  d  waves .  The  outgoing  particle  in  the  reaction 
channel  must  be  a  neutron,  as  mentioned  previously,  in  order  to  obtain 
a  sharp  cusp. 

In  conclusion,  we  remark  that  the  expressions  for  "energy- 
averaged"  Wigner  cusps  depend  only  upon  statistical  assumptions,  and 
not  upon  a  particular  model,  for  example,  the  optical  model;  rather, 
these  expressions  simply  give  the  change  in  cross  sections  expected  at 

reaction  thresholds,  due  to  unitarity  of  the  scattering  matrix.  Meyerhof 

6 
showed,   however,  that  these  statistical  assumptions  are  consistent  with 

the  optical  model  assumption  that  the  scattering  matrix  element  consists 

22 
of  a  smooth  plus  fluctuating  part. 

Also,  we  note  that  "energy-averaged"  Wigner  cusp  expressions 

predict  no  below- threshold  effect.   Consequently,  in  regions  of  mass 

number  and  energy  where  "energy -averaged"  theory  applies,  it  is  not 

necessary  to  specify  the  dominant  outgoing  wave  in  the  reaction  channel, 

in  contrast  to  pure  Wigner  cusp  theory  (see  the  discussion  just  before 

Eq.  (7))«   In  particular,  outgoing  j>  wave  dominance  of  the  reaction  in 

90 
the  Zr   experiment,  discussed  earlier,  causes  no  difficulty  in  analysis. 

Determination  of  the  dominant  incident  partial  wave  in  the  reaction  and 

Eq.  (26)  suffice. 


CHAPTER  III 
EXPERIMENT 

Method 

Neutron  Differential  Elastic  Cross 
Section  Measurements 

Neutrons  were  produced  by  the  T(p,n)  reaction,  using  tritium 
loaded  zirconium  and  titanium  targets  for  the  proton  beam  from  a  High 
Voltage  Engineering  Company  3-Mev  Van  de  Graaff  generator  stabilized 
to  about  1  kev.  Various  tritium  targets,  described  in  Table  I,  of  from 
10  to  70  kev  thickness  at  1.02  Mev  proton  energy  (as  determined  by  the 
rise  method  at  threshold)  were  used.  An  eccentric  circular  motion  was 
imposed  upon  the  tritium  target .  Consequently,  the  proton  beam  struck 
a  ring-shaped  area  on  the  target,  which  minimized  cooling  difficulties 
and  also  effectively  averaged  over  nonuniform  tritium  deposition  in 
the  target,  thereby  producing  a  time  integrated  neutron  flux  very 
nearly  proportional  to  accumulated  proton  charge . 

The  neutron  beam  was  collimated  by  a  rectangular  conical  open- 
ing in  LiC0_-loaded  paraffin  source  shielding  of  characteristic  thick- 
ness 18  in.   The  geometry  and  shielding  of  the  experiment  is  shown  in 
Fig.  1.   The  collimation  angle  was  such  that  the  energy  spread  in  the 
neutron  beam  due  to  T(p,n)  kinematics  was  approximately  k   kev.  A 
standard  long  counter   placed  at  0°  with  respect  to  the  beam  (see  Fig. 
l)  was  used  to  monitor  the  neutron  yield. 
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Fig.  1.  Experimental  Arrangement- -Plan  Views   1.  Tritium 
loaded  target;  2.  Air-water  spray  coolant;  3«  Eccentric  wobble  for 
additional  cooling  effect;  k.   Paraffin  shielding  (loaded  with 
LiCO  );  5-  Paraffin  (LiCO  loaded)  rectangular  conical  collimator; 
6.  Scatterer  in  "in  Beam"  position;  "J.   Scatterer  in  "Out"  position; 
8.  Anthracene  or  stilbene  crystal  (a  right  circular  cylinder,  2" 
dia  x  l"  long);  9*  6810A  RCA  photomultiplier  tube;  10.  Preampli- 
fier (linear  and  nonlinear  gamma  suppression  outputs);  11.  Long 
Counter  (BF  )  in  monitor  position . 
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The  neutron  beam  was  allowed  to  strike  the  scatterer,  which  was 
made  slightly  smaller  than  the  aperture  of  the  neutron  team,  at  about 
10  in.  from  the  collimator  end.   The  lithium  scatterer  was  a  right  cir- 
cular cylinder  oriented  vertically;  the  remaining  scatterers  used  were 
in  the  shape  of  a  slab,  oriented  with  their  normals  parallel  to  the 
floor  and  at  40°  to  the  incident  beam  direction.   The  scatterers  used 
are  described  in  detail  in  Table  II.  Scatterer  thicknesses  were  be- 
tween l/5  and  5/8  neutron  mean  free  paths  at  the  energies  of  experi- 
mental interest.   These  rather  thick  scatterers  were  necessary  to  obtain 
sufficient  counting  rates .   Scattered  neutrons  were  detected  at  two 
angles  by  organic  scintillation  crystals  (2  in.  dia  x  1  in.  long  anthra- 
cene and  2  in.  dia  x  3 A  in.  long  stilbene)  placed  one  on  each  of  the 
transmission  and  reflection  sides  of  the  scatterer  at  distances  varying 
from  6  to  10  in.  and  with  the  crystals  oriented  as  shown  in  Fig.  1. 
The  scatterer  to  crystal  distance  varied  with  detection  angle  and  was 
set  at  the  maximum  consistent  with  a  neutron  background/effect  ratio 
less  than  two.  A  saturated  space-charge  pulse-shape  discrimination 
output  of  the  detection  preamplifiers  was  used  to  suppress  gamma-ray 
background.   Typical  neutron  pulse  spectra,  given  in  Fig.  2,  illustrate 
that  we  also  used  an  integral  discriminator  to  bias  out  inelastic  neu- 
trons . 

To  facilitate  neutron  background  subtraction,  the  scatterers 
were  suspended  on  thin  wires  from  an  overhead  crane  so  that  the  scat- 
terer could  be  alternately  removed  and  replaced  in  the  beam  at  each 
neutron  energy  (see  Fig.  l) . 
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TABLE  I 
Description  of  Tritium- Loaded  Targets 


Index 

, 

Made  by 

Yield 

Thickness 

neutrons/sterad 
microcoulomb 
2.50  Mev 

(at  T( 
old 

p,n)  thresh - 
1.019  Mev) 
kev 

1 

HAR 

1.5  x  io6 

22 

2 

HAR 

O.37  x  10 

7 

3 

HAR 

k.Q    x  106 

60 

k 

OR 

0.33  x  10 

11 

5 

OR 

6 

1.1  X  10 

11 

BAR  =  Harwell,  England  (Tritium- loaded  titanium) 
OR  =  Oak  Ridge  Natl.  Lab.,  Tenn.  (Tritium- loaded  zirconium) 


TABLE 

II 

Description  of  Scatterers 

Scatterers  Used  in  i 

la  Measurement 

s 

Material 

Dimensions 
(inches) 

Chemical 
Purity 

Angle  between 

slab  and  beam 

(degrees) 

Cerium  a. 

k.O    x  3.0  x  O.63 

99*8 

50 

b. 

3.25  x  2.75  x  0.19 

99-8 

50 

Barium 

k.o    x  3.0  x  1.19 

97-5 

^5 

Zirconium 

k.o    x  1.5  x  0,56 

98 

50 

Iron 

k.o    x  2.5  x  0.52 

98 

50 

Lithium 

3.05  dia  X  k.k   high 

~  100 

Rt.  cir„  cylinder 

Samples  Used  in  a  (Transmission) 

Measurements 

Cerium 
Barium 
Iron 

1.0  X  1.25  x  0.935  long 
1.0  dia.  X  I.87  long 
1.0  X  1.13  X  0.97  long 

The  "long"  dimen- 
sion was  parallel 
to  the  beam  axis. 

Fig.  2.  Pulse  spectra  from  detection  equipment  at  neutron 
energy  of  0.5  Mev.  a.  Gamma-Ray-Suppression  Spectra  (nonlinear  pulses): 
1.  With  the  detector  (2"  dia  X  3 A"  long  stilbene  crystal)  in  the  neu- 
tron beam  and  with  the  pulse  height  analyzer  (PHA)  recording  all  preampli- 
fier nonlinear  output  pulses.  2.  With  the  proton  "beam  blocked  by  a  beam 

22 

stopper  from  hitting  the  tritium  target,  with  a  Na   source  taped  on  the 

crystal,  and  with  the  PHA  recording  all  preamplifier  nonlinear  output 

pulses.   3-  With  the  proton  beam  on  the  tritium  target,  with  the  detector 

22 
crystal  in  the  neutron  beam,  and  with  a  Na   source  on  the  crystal.   Two 

integral  discriminator  biases  were  set,  one  on  each  of  the  linear  and 
nonlinear  spectra,  and  the  PHA  had  been  adjusted  to  accept  nonlinear 
pulses  which  were  associated  with  a  time  coincidence  of  discriminated 
linear  and  nonlinear  pulses,  b.  Linear  Spectra  (linear  preamplifier 
output  pulses  associated  with  a  time  coincidence  of  discriminated  linear 
and  nonlinear  pulses);  A.  With  the  detector  crystal  in  a  typical  experi- 
mental position  (90°  lab)  and  the  scatterer  in  the  experimental  position 
in  the  neutron  beam.   B.  With  the  detector  in  the  same  experimental  posi- 
tion as  in  A,  and  with  the  scatterer  retracted  to  its  "Out"  position 
(see  Fig.  l) .  C.  (A  -  B)  =  elastically  scattered  neutron  pulse  spectrum. 
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Neutron  Total  Cross  Section  Measurements 

The  same  neutron  "beam,  source  shielding,  collimation,  and  gamma 
suppression  described  above  were  used  in  measuring  total  cross  sections. 

Neutrons  were  detected  in  this  case  by  a  stilbene  crystal  with 
dimensions  l/2  in.  dia  X  l/2  in.  long  mounted  at  0°  to  the  incident 
beam.   Transmission  samples  were  right  parallelepipeds  whose  area  normal 
to  the  beam  (*  1  sq.  in.)  was  sufficient  to  shield  the  crystal  from 
direct  neutrons.   (Neutrons  were  produced  over  a  proton  beam  spot  area 
of  about  l/k   in.  X  l/k   in.)  The  samples  were  10  in.  from  the  collima- 
tion (28  in.  from  the  neutron  source)  and  20  in.  from  the  crystal,, 

Sample  thicknesses  were  chosen  to  allow  approximately  60  percent  neutron 

2k 
transmission  at  energies  of  interest  in  the  experiment.    The  sample 

and  a  9  in.  long  1  in.  dia.  tungsten  shield  were  mounted  on  thin  sup- 
porting rods  above  the  ends  of  two  arms  of  a  rotatable  horizontal  cross. 
A  third  arm  had  only  a  thin  supporting  rod  attached.   Thus  it  was  pos- 
sible to  count  transmitted,  background,  and  direct  neutrons  in  sequence 
at  each  neutron  energy  by  rotating  the  cross. 

Apparatus 

It  remains  to  describe  the  electronics  systems  necessary  to  re- 
duce crystal  scintillation  events  to  (elastically  scattered)  neutron 
counts.  A  detailed  discussion  of  our  gamma  suppression  circuit  is  con- 
tained in  Appendix  D.   The  block  diagram,  Fig„  J»,   will  be  useful  in 
following  the  present  discussion. 
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Fig.  3-   Block  diagram  of  electronics,  shoving  the  function  of 
each  major  component . 
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For  each  of  two  detection  systems  an  RCA  6810A  photomultiplier 
and  preamplifier  (Figo  22)  produced  linear  and  nonlinear  (gamma  sup- 
pression) outputs  in  response  to  a  crystal  scintillation .  These  two 
outputs  fed  into  a  double  delay  line  and  RC  coupled  linear  amplifier, 
respectively. 

Each  of  the  amplifier  outputs  drove  an  adjustable  bias  circuit 
(integral  discriminator)  which  generated  a  standard  square  20  volt 
output  for  each  input  pulse  that  exceeded  the  bias  level o  After  time 
synchronization  was  established  by  the  use  of  a  suitable  delay  cable, 
these  pulses  were  allowed  to  enter  a  standard  coincidence  unit.  A  co- 
incidence output  pulse  was  obtained  for  each  simultaneous  occurrence 
of  input  pulses.   Thus,  with  properly  adjusted  bias  levels  on  each  of 
the  linear  and  nonlinear  pulses,  a  coincidence  output  pulse  signalled 
the  detection  of  an  elastically  scattered  (linear  pulse  above  linear 
bias)  neutron  (nonlinear  pulse  above  nonlinear  bias).,  These  pulses 
were  the  experimental  counts „ 

A  256-channel  pulse-height  analyzer  (PRA)  was  a  valuable  tool 
in  setting  up  the  gamma  suppression  circuit,  the  amplifier  gains,  and 
the  discriminator  biases.  We  also  used  the  PHA  as  a  monitor  on  the 
overall  electronics  performance  during  a  run;  these  monitor  spectra  had 
a  specific  role  in  data  analysis  for  some  of  the  experiments  (to  be  dis- 
cussed). Typical  pulse  spectra  are  shown  in  Fig.  2.   The  coincidence 
output  pulses  also  served  to  gate  the  PHA  when  spectra  such  as  those 

shown  in  Fig.  2  were  taken. 

23 
A  standard  Hansen  and  McKibben  long  counter   was  employed  to 
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monitor  the  neutron  yield.   Direct  neutrons  were  moderated  by  a  paraffin 

10  3 

cylinder  and  some  were  captured  by  boron   in  a  one-in.  dia.  BF  detec- 
tion tube.   The  subsequently  emitted  alpha  particle  initiated  a  pulse 

which  was  amplified  and  counted  as  a  neutron  event.  The  efficiency  of 

-3 
the  counter  was  approximately  2  X  10   and  was  fairly  constant  for  a 

range  of  neutron  energies  from  0.5  Mev  to  3  Mev.  A  l/8  in.  cadmium  and 

3  in.  paraffin  shield  on  the  sides  and  back  of  the  counter  reduced 

floor  scattered  neutrons  to  less  than  15  percent  of  direct  neutrons „ 

Although  absolute  neutron  yields  of  our  tritium  targets  were  measured 

using  the  long  counter  (see  Table  i),  these  yields  were  not  required 

in  the  data  analysis . 

Procedure 

Differential  Elastic  Cross  Section 

Since  the  experimental  interest  centered  on  effects  near  the 
energies  of  first  excited  states,  initial  adjustments  of  the  apparatus 
were  made  at  these  energies .  For  each  experiment  the  nonlinear  preamp- 
lifier controls  (see  Appendix  D)  were  adjusted  for  optimum  gamma-ray 
suppression,  while  the  crystal  was  located  in  the  scatterer  "in  beam" 
position,  later  called  the  calibration  position o"  Care  was  taken  to 
insure  that  the  crystal  plane  was  parallel  to  the  incident  beam  axis. 
This  aspect  corresponded  to  the  aspect  of  the  crystal  toward  scattered 
neutrons  when  it  was  in  experimental  position  (see  Fig.  l)o  After 
these  adjustments  were  completed,  the  neutron  yield  and  overall  effi- 
ciency of  the  crystal/electronics  system  were  calibrated  by  counting 
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direct  target  neutrons  at  several  neutron  energies  above  and  "below 
thresholds   This  was  done  for  each  of  two  crystal  detector  systems. 

After  the  crystals  were  moved  to  the  experimental  position 
(Fig,  l),  the  same  energy  range  was  covered  in  energy  increments  of 
approximately  one-half  incident  neutron -energy  spread.  At  each  energy, 
counts  were  recorded  both  with  the  scatterer  in  and  with  the  scatterer 
out  of  the  beam,  thus  providing  a  basis  for  subtraction  of  neutron 
background.   On  a  succeeding  run  down  in  energy  the  same  increments 
were  used,  but  the  points  interleaved  those  from  the  run  up;  thus,  we 
checked  data  reproducibility. 

All  counts  were  normalized  to  a  fixed  amount  of  accumulated 
proton  charge.  When  using  30-50  microamperes  proton  current,  we  found 
that  the  necessary  counting  time  was  typically  10  minutes  for  each 
energy  point. 

Long  counter  monitor  neutron,  counts  were  also  recorded  at  each 
point.  Monitor  neutron  counts  with  the  scatterer  out  of  the  beam 
served  both  as  a  check  on  the  condition  of  the  tritium  target,  which 
was  subjected  at  times  to  as  much  as  70  watts  of  heat  flow,  and  as  the 
basis  for  a  small  correction  of  the  data  analysis  for  target  yield 
variation. 

Total  Cross  Sections 

As  stated  above,  the  gamma  suppression  circuit  was  adjusted  at 
a  neutron  energy  corresponding  to  an  excited  state  of  interest.  After 
carefully  aligning  the  sample  and  shield  between  the  proton  beam  spot 
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and  crystal  detector,  we  took  a  sequence  of  counts  at  each  energy  point, 
By  rotating  the  supporting  horizontal  cross  on  which  sample  and  shield 
were  mounted  we  were  able  to  place  first  the  sample,  then  the  shield, 
and  last  the  duplicate  sample  support  (for  a  count  of  direct  neutrons) 
into  the  beam.  After  each  such  set  of  counts  the  neutron  beam  energy 
was  increased  by  approximately  one-half  incident  beam-energy  spread 
until  a  desired  energy  range  was  covered  below  and  above  threshold.  A 
subsequent  run  down  in  energy,  with  interleaving  points,  verified  re- 
producibility . 


CHAPTER  IV 
DATA  REDUCTION  AND  ANALYSIS 

Data  Normalization 

The  basic  normalization  used  in  the  experiment  was  to  a  fixed 
accumulated  proton  charge,  which  implies  an  assumption  that  the  result- 
ing time-integrated  neutron  flux  was  also  a  fixed  quantity,  except  for 
the  smooth  change  in  the  T(p,n)  reaction  cross  section  with  energy, 
and  except  for  a  ±1  percent  uncertainty  introduced  by  our  proton  cur- 
rent integrator.  As  mentioned  before,  the  eccentric  circular  motion 
of  the  tritium  target  minimized  any  imperfection  in  this  assumption, 
by  averaging  over  nonuniformities  of  tritium  deposition.  However, 
small  variations  in  integrated  neutron  flux  from  the  true  smooth  shape 
of  the  T(p,n)  reaction  cross  section  might  have  been  present  when  the 
neutron  energy  was  changed,  because  of  changes  in  proton  beam  intensity 
distribution  which  occur  when  one  sets  the  accelerator  to  a  new  energy. 
If  corrections  were  not  made  this  could  have  produced  scatter  in  the 
data,  as  plotted  versus  energy,  in  excess  of  statistical  fluctuations. 
For  this  reason,  monitor  neutron  counts,  which  had  very  small  statis- 
tical errors,  were  used  to  renormalize  the  differential  cross  section 
effect  counts  to  a  smooth  curve  drawn  through  monitor  neutron  data; 
this  curve  was  assumed  to  have  the  true  smooth  shape  of  the  neutron 
yield  versus  neutron  energy.   The  normalization  was  then  a  percentage 
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correction  to  each  data  point  depending  upon  the  percent  deviation  of 
the  monitor  counts  from  the  smooth  curve.   This  correction  was  usually 
less  than  one  percent  and  never  more  than  three  percent.   It  was  gen- 
erally impossible  to  distinguish  by  eye  the  difference  between  corrected 
and  uncorrected  data,  but  in  seeking  cusp  effects,  we  wished  to  take 
all  precautions  necessary  to  establish  the  relative  magnitude  of  the 
data  points  plotted  versus  neutron  energy. 

Data  Averaging 

The  original  data  for  experiments  on  zirconium,  cerium  and 
barium  showed  reproducible  fluctuations  in  addition  to  a  well  defined 
drop  in  magnitude  over  a  broad  region  above  reaction  thresholds .   It 
turned  out  that  energy  averaging  the  original  data  did  not  mask  the 
essential  structure  or  features.  Figure  k   illustrates  this  point, 
where  we  have  used  differential  elastic  cross  section  data  from  a 
cerium  experiment.   Therefore,  we  were  able  to  improve  the  statistical 
errors  without  significant  loss  of  structural  detail  and  also  to  estab- 
lish a  standard  energy  averaging  interval  (A  E  =50  kev)  for  presenting 
measurements  at  many  angles . 

Insofar  as  we  were  concerned  with  pure  Wigner  cusp  effects,  we 
could  not  use  data  averaging  procedures  in  the  experiments  on  lithium 
and  iron.  Here,  the  magnitude  of  the  structure  in  the  cross  sections 
turned  out  to  exceed  that  expected  of  cusp  effects;  in  the  experiments 
on  iron,  the  structure  was  particularly  marked.  A  successful  search 
for  a  cusp  effect  would  depend  upon  distinguishing  this  basic  structure 
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Fig.  k.      Neutron  differential  elastic  (90°)  cross  section  of 
natural  cerium  for  three  widths  of  energy  averaging  interval, 
(a)  Original  data;  A  E  =  7  kev;  statistical  errors  ±5$.   (b)  Same  data 
averaged  over  A  E  =  25  kev;  statistical  errors  2.5$.   (c)  Result  of  a 
larger  averaging  interval,  A  E  =  50  kev;  statistical  errors  approxi- 
mately 1$   The  averaging  in  (c)  was  selected  as  the  standard  for  cusp 
analysis  of  cerium.   The  loss  of  structural  detail  upon  successive 
averaging  is  apparent. 
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in  the  cross  sections  from  a  true  cusp  effect.  We  illustrate  with 
Fig.  5,  which  gives  iron  total  cross  section  measurements,  and  which 
shows  that  use  of  25  kev  incident  neutron  energy  resolution  washed 
out  significant  structure  elucidated  by  experiments  with  7  kev  spread 
in  incident  neutron  energy. 

In  addition  to  washing  out  structure,  a  poor-resolution  inci- 
dent beam  also  has  the  well  known  effect  of  distorting  sharp  features 
in  cross  sections.  We  observed  this  in  differential  cross  section 
measurements;  a  comparison  of  Fig.  11  and  Fig.  Ik   shows  a  feature  shift 
of  about  +10  kev  as  the  result  of  using  30-kev  compared  to  7-kev  energy 
resolution.   This  amount  of  feature  shift  would  make  a  cusp  analysis 
impossible. 

In  reporting  some  of  the  lithium  differential  cross  sections 
we  have  energy  averaged  the  data  for  clarity  in  presentation  (see  Fig. 
10) .   It  is  clear  from  the  foregoing  discussion  that  in  so  doing  we 
conceded  that  there  was  no  cusp  effect  in  the  data  and  that  remaining 
structure  must  be  attributed  strictly  to  the  elastic  scattering  cross 
section  of  the  lithium  nucleus. 

Data  Analysis 

Differential  Cross  Sections 

We  give  our  calculation  of  the  measured  differential  elastic 
cross  sections  using  (l)  the  basic  data,  which  consisted  of  the  counts 
C  =  CT  ;{scatterer  in)  -  C   (scatterer  out),  (2)  geometrical  factors, 


kk 


Fig.  5-  The  total  neutron  cross  section  of  natural  iron,  using 
neutron  beams  of  two  energy  spreads:   (a)  A  E  =  25  kev;  (b)  A  E  =  7 
kev.   The  size  of  the  symbol  indicates  the  size  of  the  statistical 
errors . 
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and  (3)  an  overall  efficiency  measured  in  a  calibration  run.   In  order 
to  make  a  detailed  comparison  of  gross  structure,  cusp  effects  in 
cross  sections,  or  angular  distributions  with  theory,  or  with  other 
experimental  work,  additional  corrections  for  multiple  scattering  and 

finite  scatterer  effects  had  to  be  applied.  These  corrections  are  dis- 

) 

cussed  in  Appendix  E  and  tabulated  in  Table  III. 

Measured  cross  sections  will  be  given  in  the  center  of  mass 
(CM)  frame  and  will  be  expressed  as  functions  of  mean  polar  CM  angle 
(0  (CM))  and  mean  CM  incident  neutron  energy  (E  (CM)).   Our  reason 
for  presenting  CM  data  is  that  theoretical  formulas  for  cusp  phenomena 
and  threshold  energies  are  more  conveniently  expressed  in  CM  terms . 
The  quantities  da(CM),  0(CM),  and  E(CM)  were  computed  from  the  cor- 
responding quantities  expressed  in  the  laboratory  frame,  using  well 

25 

known  transformation  formulas.  In  the  present  experiment  the  differ- 
ences between  lab  (laboratory)  and  CM  quantities  were  negligible  except 
in  the  cases  of  iron  and  lithium. 

E(lab)  was  determined  by  subtracting  one-half  the  incident 
neutron  energy  spread  from  the  maximum  lab  energy  of  the  neutron  beam. 
Incident  neutron  energy  spread  was  computed,  assuming  statistical  in- 
dependence of  the  two  contributions,  from  (l)  the  geometrical  spread, 
fixed  by  the  neutron  collimation  angle  and  T(p,n)  kinematics  to  be 
k   kev,  at  1.5-Mev  neutron  energy;  and  (2)  the  proton  energy  spread 
caused  by  the  thickness  of  the  metal  target,  containing  tritium.  The 

latter  spread  was  computed  from  the  target  thickness  at  T(p,n)  thresh - 

23 
old  "  [l.019-Mev  proton  energy  (lab)]  as  determined  from  the  slope  of 


hi 


the  rise  at  threshold  and  conversion  factors  relating  proton  stopping 

power  of  the  target  at  threshold  to  the  stopping  power  at  the  energy 

25 
of  interest.    Neutron  energy  spread  is  directly  related  to  proton 

/     N  26 

energy  spread  by  T(p,n;  kinematics „ 

Maximum  incident  neutron  lab  energy  present  in  the  incident 
beam  was  determined  from  maximum  proton  lab  energy  by  the  same  kine- 
matical  relations,  where  the  proton  energy  was  calibrated  in  two  ways. 
First,  by  using  the  accurately  known  T(p,n)  threshold  energy  we  ef- 
fectively calibrated  the  Van  de  Graaf  machine  analyzing  magnet  at 
1.019  Mev.   The  magnet  field  was  measured  with  a  proton  magnetic  res- 
onance probe,  between  the  pole  faces,  to  obtain  with  high  precision 
an  RF  frequency  which  could  be  associated  with  1.019  Mev  proton  energy. 
We  used  the  intersection  of  the  straight  line  extrapolation  of  the 
lead  slope  of  the  T(p,n)  yield  versus  frequency  curve  and  the  frequency 
axis  as  indicating  the  true  threshold.   Since  proton  resonance  fre- 
quency and  the  magnetic  field  (B)  are  linearly  related  (f  =  u  °  B/h, 
where  f  is  frequency,  u  is  the  proton  magnetic  moment,  and  h  is 
Planck's  constant),  and  the  analyzing  magnet  is  a  momentum  selector, 

then  maximum  proton  lab  energy  is  B   (lab)  =  1.0 19  f  /fmu  Mev  where 

p  In 

2 
f     is  the  previously  described  threshold  proton  resonance  frequency. 
In 

12  s 

Second,  because  C    has  a  sharp  resonance  at  2 .07o-Mev  lab 

neutron  energy,  a  long  counter  of  the  type  used  in  this  experiment  has 

an  anomaly  in  its  efficiency  versus  energy  curve  at  2.076  Mev  due  to 

12  23 

the  increased  cross  section  of  C    in  the  paraffin  moderator.    Thus, 

a  dip  occurs  in  monitor  neutron  counts  at  this  neutron  energy.   The 
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proton  resonance  frequency  at  this  dip  established  a  calibration  point 
for  the  mean  neutron  energy  which  served  as  a  check  on  calibration  by 
the  first  method.   The  two  calibration  points  conveniently  straddled 
the  proton  energy  region  involved  in  the  present  experiments. 

Initially,  the  "effect  counts"  C  were  normalized  as  described 
in  the  first  section  of  this  chapter,  giving  normalized  effect  counts, 
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For  a  given  scattering  sample  and  incident  neutron  lab  energy 


E  the  relation  between  C„  and  the  differential  cross  section  is: 

L-If  [(I-L)jH]  fjE  -  Ej  fl  fl 
N      x  v     R"   x   Dv     R   x  s 


X  [1  -  exp(-  naT  d)]  da(0,  E)/oT  (27) 

where  Y  is  the  neutron  yield  of  the  target  per  unit  solid  angle; 

f  [(E  -  E  ),  R  ]   is  the  intrinsic  crystal  efficiency  at  the  lab  energy 
x       K    x 

of  the  scattered  neutron  and  at  the  crystal-to-scatterer  distance,   R  ; 

E   is  the  lab  recoil  energy  of  the  residual  nucleus;  f^E  -  E  )  is 
R  D      R 

a  discriminator  correction  to  be  discussed:   ii   and  0,       are  the  solid 

'   x       s 

angles  subtended  by  the  crystal  to  the  scatterer  center  and  by  the 
scatterer  to  the  proton  beam  spot,  respectively;   nd  is  the  number  of 
scattering  nuclei  per  unit  area  in  the  direction  of  the  incident  beam; 
aT  is  the  total  cross  section  of  the  scatterer  material;  and  dcr(0,  E) 
is  the  differential  elastic  cross  section  per  unit  solid  angle.   Since 
dcr(0,  E)  has  been  assumed  constant  [otherwise  Eq.  (27)  would  contain 
integrals  over  solid  angle  of  the  finite  detector  and  energy  spread  of 
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the  beam]  over  the  angular  region  described  by  £2   and  the  energy 
region  included  by  the  energy  spread  of  the  incident  neutrons,  the 
symbols  6     and  E  refer  to  mean  polar  angle  and  mean  neutron  beam 
energy  respectively  [i.e.,  6   (lab)  and  E(lab)]„ 

The  corresponding  formula  for  "calibration  counts"  is 

C   =  Y  fl  f  (E,R  )  fJE),  (28) 

c      c  x    c   D 

where  Q       is  the  solid  angle  subtended  by  the  crystal  to  the  proton 
beam  spot  and  R   is  the  associated  scatterer  to  tritium  target  dis- 
tance. 

Taking  the  ratio  Eq.  (27)/Eq.  (28)  yields 

f  (E,R  )  fD(E)  ft 

d(j(0'E)  =  f  [(E  -  Ej,  R  ]  f  (E  -  E_)  ?fV 
x       R    x   D      R    x  s 

X  oT[l  -  exp(-  n  aT  d)]"1  K(CN/Cc)  (29) 

where  K  is  a  constant  normalizing  the  effect  and  calibration  counts 
to  the  same  accumulated  proton  charge. 

The  factors  f   are  necessary  corrections  which  arise  because 
the  discriminator  bias  on  detected  linear  pulses  was  fixed .  As  a  re- 
sult, the  counted  fraction  of  the  total  detected  elastic  neutron  spec- 
trum is  a  function  of  energy.   This  is  shown  in  detail  by  Fig.  6  and 
is  explained  further  in  the  figure  caption.   In  the  detection  of  scat- 
tered neutrons,  the  factor  f^  becomes  fp,(E  -  E  )  because  the  lab 
'  D  D      R 

energy  of  the  scattered  neutrons  rather  than  that  of  the  incident  neu- 
trons determines  the  discriminator  correction.   In  all  except  the 
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Fig.  6.   Linear  pulse  spectra  for  elastically  scattered  neutrons 
at  three  energies,  shoving  the  energy  dependence  of  detection  efficiency 
due  to  fixed  integral  discrimination  bias.   The  spectra  were  obtained 
with  the  detector  at  125°  (lab)  to  a  lithium  scatterer  and  are  the  re- 
sult of  subtracting  the  spectrum  with  scatterer  out  of  the  neutron  beam 
from  the  spectrum  with  scatterer  in  the  beam.  This  figure  only  illus- 
trates an  experimental  point  in  the  data  analysis.  Actual  runs  were 
made  with  the  linear  bias  set  at  a  fixed  voltage,  corresponding  to  ap- 
proximately l/2  maximum  pulse  height  for  a  neutron  energy  equal  to  the 
mean  energy  of  a  complete  run.   In  this  example,  the  bias  would  be  set 
at  about  pulse  height  70,  but  the  run  could  not  extend  as  low  as  0.28- 
Mev  neutron  energy. 
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experiment  on  LI  and  Fe  f^E  -  E  )  *»  f-^E)  because  the  lab  recoil 

D  K      1) 

energy  of  the  residual  nucleus  was  negligible.  Thus  the  ratio 

f_(E)/f_(E  -  E  )  in  Eq.  (29)  was  *»1.  For  the  lithium  and  iron  experi- 
D     D      R 

ments  we  determined  this  ratio  by  comparison  of  pulse  height  analyzer 
spectra  of  linear  pulses  taken  during  "calibration"  and  "effect  runs. 
We  simplified  the  ratio  f  (E,  R  )/f  [ (E  -  E  ),  R  ],     which 

.X.  (_■      A\  X\        A. 

is  a  ratio  of  intrinsic  crystal  efficiencies,  in  a  reasonable  way,  by 
assuming  when  necessary  (in  lithium  and  iron  experiments)  that  the 
energy  dependence  of  f   was  separable  from  the  distance  dependence, 

A 

R  .  Recalling  the  calibration  equation,  Eq.  (28),  we  see  that 

Cc/[Y(E)  fic  fD(E)]  =   fx(E,RQ)  (30) 

Thus,  if  we  divide  the  plot  of  C   versus  E  by  the  incident  flux, 
Y(E),   obtained  from  monitor  neutron  counts,  and  by  f  (E),   obtained 
from  linear  pulse  spectra,  we  see  that  the  energy  dependence  of  f 

A 

is  determined.   The  same  energy  dependence,  by  assumption,  applied  to 

f  (E  -  E  ,R  ),   hence  by  translation  of  the  curve  f   along  the  energy 

axis  by  an  amount  E   (calculated  from  kinematical  relations   )  we 

R 

obtained  the  ratio  f(E  -  E  ,R  )/f  (E,  R  )  insofar  as  the  (separated) 

R   X    X      C 

energy  dependence  is  concerned,    i.e., 

f  (E,    R   )/f   [(E  -  Ej,    R   ]    =     €),(E)   f   (R   )/f   (R   ),  (31) 

X  Q.     '     X  R"        X  V      '        XV     Q."     X       X"  ' 

where  €. (E)  is  the  just  determined  energy  dependence  of  the  f  ratio, 
It  will  be  useful  to  write 

%/ax  %  ,  Sl(H  2/Aj,  (32) 
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where  A   is  the  scatterer  area  projected  perpendicularly  to  the  beam 
s 

axis  and  e,   is  a  small  correction  for  the  degree  to  which  the  three 
solid  angles  are  not  exactly  expressible  as  projected  area/distance 
squared. 

Substituting  (31)  and  (32)  in  Eq.  (29)  we  obtain 

f  ( R  )  f  (E ) 

da(9,E)  =  %<l)  j^j   ei(Rx2/As)  fp(g  .  Er) 

X  aT(l  -  exp(-  n  a^   d))'1  K  CN/Cc  (33) 

The  quantity  e,  f  (R  )/f  (R  )  =  e9(R  )  represents  the  inverse 
of  the  effective  crystal  efficiency  as  a  function  of  distance,  R  , 
normalized  to  its  value  at  R  =  R  .  If  we  make  the  reasonable  assump- 

.A.        *— 

tion  that  this  quantity  is  very  close  to  unity  for  R   equal  to  and 

greater  than  R   (R  was  28  in.  in  our  experiment)  then  neutron  counts 
c    c 

versus  distance  can  be  compared  to  inverse  square  law  behavior  (normal- 
ized at  large  distances)  to  obtain  ep(R  ).  This  measurement  was  made 
using  a  PoBe  neutron  source  at  several  distances  (R  )  from  the 
crystal,  to  obtain  relative  counts  in  unit  time  for,  each  position. 
The  results  are  shown  in  Fig.  7«  Large  corrections  would  be  necessary 
for  distances  less  than  5  in.  In  the  present  experiments  the  factor 
€p(R  )  was  typically  1.10  to  1.20  (R  varied  from  5.5  to  9  in.  de- 
pending upon  the  polar  angle  6     of  the  cross  section  measurement). 

In  view  of  the  foregoing  reductions  we  may  write  the  most  use- 
ful form  of  the  expression  for  cross  section: 


5h 


_2 
Fig,  7-  The  deviation  from  (distance)  "  efficiency  response  of 

a  2"  dia  X  1  long  stilbene  crystal,  oriented  with  a  diameter  parallel 

to  the  distance  coordinate.   Curve  A:   (Distance)   fitted  to  the  data 

with  largest  abscissa.  Curve  B:  Smooth  curve  through  the  experimental 

points.  Curve  C:  Calculated  curve  using  graphical  analysis  and  effi- 

23 
ciency  corrections  from  Marion  and  Fowler. 
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x  aT(l  -  exp(-  n  aT  d))"1  K  CN(E)/Cc(E)  (3*0 

where  we  have  included  a  factor  €   denoting  the  conversion  of  the 
cross  section  to  the  CM  frame  -  For  all  except  the  experiments  on  iron 
and  lithium,  the  factor  e, (E)  f  (E)/f  (E  -  E  )  was  unity,  which 
simplified  data  analysis . 

Total  Cross  Sections 
The  data  analysis  for  total  cross  sections  is  described  by 
aT(E(CM))  =  e^ND)"1  ln[(CD  -  Cg)/(CT  -  Cfi)]  (35) 

where  aT(E(CM))  is  the  total  cross  section  as  a  function  of  mean 
center  of  mass  neutron  energy,  e       is  (as  before)  a  symbol  of  the 
conversion  of  the  energy  dependence  from  lab  to  CM  frame,   ND  is  the 
number  of  nuclei  per  unit  area  in  the  sample,   In  means  natural  log- 
arithm, C   is  the  number  of  direct  neutron  counts,  C_.  is  the  number 
of  background  counts,  C„  is  the  number  of  transmitted  counts,  and 
where  C_.,  C  ,  and  Cm  are  all  functions  of  neutron  energy  and  are 

Do  1 

all  normalized  to  a  fixed  accumulated  proton  charge.  No  monitor  neu- 
tron normalization  was  applied  because  only  count  ratios  are  involved 
with  the  machine  energy  at  one  setting. 
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Error  Analysis 

Differential  Elastic  Cross  Sections 

Contributions  to  the  uncertainty  in  the  absolute  magnitude  of 
the  differential  elastic  cross  sections  were  present  in  each  factor  of 
Eq..  (3*0>  reproduced  below; 

dWeCM'fN(CM)>  -  ^xX)  €2<V  €3  €U(E)  (fD(E)/fD(E  "  V} 

X  aT(l  -  exp(-  n  a  d))"1  K  CN/Cc  (3*0 

Statistical  errors  entered  the  results  through  the  factor 

C„/C  .  Those  in  C   were  negligible .   Those  which  were  in  C„  are 
W    c  c        &  &  N 

separately  described  in  the  figure  captions  for  each  set  of  experimental 
data.   In  general,  the  statistical  errors  were  rather  small  (from  1  to 
3  percent);  hence,  they  had  little  effect  on  the  precision  of  the  abso- 
lute differential  cross  section  measurements. 

Experimental  uncertainties  in  the  factors  of  Eq.  (3^)  are 
discussed  below: 

1.  Factor  C„/C  .   There  were  two  experimental  uncertainties 

in  C  .   The  detector  was  located  at  the  scatterer  "in  beam"  position 
c 

during  the  calibration  run.   Its  distance  from  the  neutron  source  was 
about  28  in.   The  precision  with  which  the  crystal  could  be  placed  at 
the  exact  scatterer  position  was  ±l/8  in.   Hence,  there  was  an  ±1 
percent  uncertainty  from  this  source.  Also,  we  include  in  C   an 
error  associated  with  the  precision  to  which  the  crystal  was  oriented 
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with  respect  to  the  incident  neutrons.  The  aspect  of  the  crystal  rela- 
tive to  the  scattered  neutrons,  when  the  crystal  was  in  detection 
position  (see  Fig.  l),  should  have  been  the  same,  because  deviations 
produce  differences  in  effective  crystal  efficiency.  We  assign  an  ad- 
ditional uncertainty  of  ±2  percent  to  C„,   as  a  measure  of  the  preci- 
sion with  which  the  crystal  aspects  were  adjusted. 

2 

2.  Factor  (R  /A  ) •   The  experimental  uncertainty  in  R  , 

X     S  X 

the  crystal  to  scatterer  distance,  was  ±1/8  in.  for  mean  R   about  7 
in.;  thus,   R   had  an  extreme  uncertainty  of  ±2.5  percent.   The  dimen- 

X 

sions  and  aspect  (with  respect  to  the  beam  axis)  of  the  scatterers  were 

not  precisely  known,  which  led  to  a  ±1  percent  uncertainty  in  A  .  A 

s 
2 
total  ±6  percent  uncertainty  was  assigned  factor  R  /A  . 

X     o 

3.  Factor  a (l  -  exp(-  n  a  d))   .   Uncertainties  in  the  value 
of  a   enter  this  factor  only  in  second  order  terms  in  n  a  d.   The 
precision  to  which  the  factor  nd  was  known  determined  the  experimental 
error.  A  ±1  percent  uncertainty  was  assigned  this  factor. 

k.     Factor  K.  Although  the  current  integrator  used  in  deter- 
mining the  normalization  to  fixed  proton  charge  is  rated  at  ±1  percent 
reproducibility,  the  mean  normalization  for  several  counting  periods 
of  the  same  accumulated  charge  was  quite  precise.   So  also  was  the  mean 
normalization  for  several  counting  periods  of  a  different  fixed  proton 
charge  (i.e.,  the  counting  period  for  calibration  points  was  much  less 
than  that  for  runs  involving  detection  of  scattered  neutrons).   The 
uncertainty  in  K  was  at  most  ±0 . 5  percent . 
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5.  Factor  e0(R  ).  This  factor  was  obtained  by  comparing  ex- 
perimental  points  with  a  theoretical  (distance)  '  behavior  for  effective 
efficiency  versus  distance  (see  Fig.  7).  Figure  7  shows  that  the  data 
did  not  precisely  fix  the  normalization  of  the  curve  at  large  distances . 
The  figure  also  shows  some  discrepancy  between  experimental  and  calcu- 
lated points  (discussed  elsewhere  in  this  report).  Considering  R  =7 
in.  typical,  we  estimated  ±3  percent  uncertainty  in  ep(R  ). 

6.  Factor  e.,.   This  factor  is  symbolic,  representing  conver- 
sions involved  in  quoting  results  in  terms  of  mean  CM  energies,  CM 
angles,  and  CM  intensities,  where  CM  is  center  of  mass.   Quantities 

M, /M   and  M  /(M  +  M^),  which  are  precise  in  the  present  usage, 
dominate  these  conversion  formulas,  the  uncertainty  in  lab  angles  and 
energies  playing  a  minor  role.  An  extreme  uncertainty  ±0.5  percent 
was  assigned  factor  e_. 

7.  Factor  e. (E)  f_(E)/f_(E  -  E_)„  In  all  except  the  lithium 

4     D     D      K 

and  iron  experiments,  this  factor  is  unity  to  good  precision.  In 
lithium  and  iron  cross  sections  it  is  a  source  of  additional  uncertainty, 
These  uncertainties,  which  are  potentially  large  for  a  single  data 
point,  are  reduced  if  one  observes  a  smooth  trend  in  the  factor  when  it 
is  determined  for  points  at  many  energies.  We  estimated  a  ±2  percent 
uncertainty  for  iron  cross  sections  and  ±5  percent  uncertainty  for 
lithium  cross  sections  from  this  factor . 

Some  of  the  experimental  errors,  given  in  1  through  7>  were 

2 
correlated;  for  example,  R    is  anti- correlated  to  e0(R  ).  For  the 

purposes  of  this  experiment,  detailed  examination  of  such  correlations 


6o 


was  not  necessary  and  the  procedure  adopted  was  to  assume  independence 
of  the  errors;  thus,  the  sum  of  the  above  percent  uncertainties  applied. 

We  concluded  that  exclusive  of  statistical  errors  and  those 
errors  associated  with  multiple  scattering  corrections  (see  Appendix 
E),  the  absolute  differential  cross  sections  of  zirconium,  barium,  and 
cerium  were  known  to  within  ±lh   percent;  of  iron  to  within  ±l6  percent; 
and  of  lithium  to  within  ±19  percent. 

Total  Cross  Sections 

Measurements  of  total  cross  sections  were  accomplished  in  a 
neutron  background  of  less  than  5  percent  of  direct  neutron  counts  which 
minimized  second  order  background  corrections  to  the  data.  The  samples 
had  small  dimensions  perpendicular  to  the  incident  beam  direction,  be- 
ing just  sufficient  to  completely  shield  the  detector  from  direct 
neutrons.   The  sample-to-detector  and  sample- to- source  distances  were 
large  compared  to  the  detection  crystal  dimensions.   Neutron  transmis- 
sion through  the  samples  was  approximately  60  percent  of  the  incident 

flux.  All  these  are  excellent  experimental  conditions  for  measurement 

2k 
of  total  cross  sections.    A  crude  calculation  indicated  that  in- 
scattered  neutron  corrections  to  the  cross  section  were  at  most  2  per- 
cent 0   The  main  constraint  on  the  amount  of  in-scattering  was  the 
small  solid  angle  the  detector  crystal  presented  to  in-scattered  neu- 
trons. An  uncertainty  of  ±3  percent  was  assigned  to  account  for  second 
order  background  scattering  effects  and  sample  in- scattering  corrections 
which  we  did  not  calculate  in  detail.   Thus,  we  are  assuming  in  effect 
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that  the  uncertainty  in  our  rough  calculation  is  larger  than  the  cor- 
rection .  A  ±1  percent  extreme  error  was  estimated  for  the  measurement 
of  ND,   the  number  of  scattering  nuclei  per  unit  area. 

Hence,  exclusive  of  statistical  errors  which  are  given  sepa- 
rately for  each  experiment,  the  total  uncertainty  in  our  absolute  total 
cross  section  measurements  was  ±k   percent. 


CHAPTER  V 
EXPERIMENTAL  RESULTS  AND  COMPARISON  WITH  THEORY 

Energy  and  Angle  Dependence  of  the  Cross  Sections 

Although  the  basic  purpose  of  the  experiment  was  to  search  for 
threshold  effects  in  neutron  cross  sections,  data  such  as  angular  dis- 
tributions and  differential  elastic  and  total  cross  sections  were 
obtained  which  are  interesting  in  their  own  right,  exclusive  of  possible 
cusp  effects.  These  data  will  be  discussed  briefly,  with  emphasis  on 
illustrations  to  clarify  the  discussions. 

The  measured  angular  distribution  of  1„6-Mev  neutrons  scattered 
elastically  from  natural  cerium  is  presented  in  Fig.  8.   The  data 
points  are  cross  plotted  from  measurements  as  a  function  of  energy  at 

separate  angles.  As  shown  in  the  figure,  an  optical  model  prediction 

17 
of  Moldauer,   modified  for  multiple  scattering  and  finite  scatterer 

effects  (see  Appendix  E)  agrees  very  well  in  shape  and,  except  at  back 

angles,  in  magnitude,  with  this  data. 

Figure  9  shows  the  measured  angular  distribution  of  600  kev 

(CM)  neutrons  elastically  scattered  from  natural  lithium.   The  data  of 

this  experiment  compares  favorably  with  that  from  a  compilation  of 

27 

other  experiments  by  Howerton,   which  is  also  shown.   This  substan- 
tiates our  estimate  that  multiple  scattering  corrections  were  small  for 
the  lithium  experiments  (see  Appendix  E) .   The  shape  of  the  distribution 
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Fig.  8.  Angular  distribution  of  neutrons  elastically  scattered 

17 
from  natural  cerium  at  E  =1.6  Mev.  Curve  1;  A  Moldauer   optical 

n 

model  prediction.  Curve  2;  Curve  1  modified  with  a  multiple  scattering 
correction.  Curve  3:  Curve  2  modified  with  calculated  experimental 
finite  scatterer  effects.   Triangles  are  experimental  points  from  the 
present  experiment. 
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Fig.  9   Angular  distribution  of  neutrons  elastically  scattered 

from  natural  lithium  at  E  =0.6  Mev„   Open  circles  are  data  from  the 

n 

present  experiment.  This  data  is  not  corrected  for  multiple  scattering. 

27 
Triangles  are  data  from  Howerton's  compilation. 
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is  interesting,  being  peaked  backwards,  presumably  as  a  result  of  inter- 

28 
ference  between  the  strong  p_  wave  resonance   at  0.25  Mev  and  an  s 

wave  background.   This  backward  peaking  means  that  one  should  look  at 

forward  angles  for  cusp  effects,  because  the  theoretical  percentage  ef- 

1/2 
feet  is  inversely  proportional  to  (da)    (see  Chapter  II). 

In  Fig.  10  we  show  some  of  the  differential  cross  section  data, 

versus  energy,  from  which  the  angular  distribution,  Fig.  9>  was  cross- 

28 
plotted.  At  the  lower  energies  the  upper  edge  of  the  p_  wave  resonance 

at  0.25  Mev  is  in  evidence.  The  wavy  structure  in  all  the  cross  sec- 
tions is  probably  not  an  experimental  effect;  it  is  quite  surprising, 
because  for  a  light  nucleus  such  as  lithium  one  might  expect  rather 
smoothly  varying  cross  sections  at  these  energies,  except  for  isolated 

resonances.  The  inelastic  cross  section  to  the  first  excited  state  of 

7  29 

Li  ,  as  measured  by  Freeman,  Lane,  and  Rose   shows  similar  structure. 

The  differential  elastic  cross  section  of  iron,  versus  energy, 

was  measured  for  four  angles  using  a  rather  thick  neutron  target  (#3> 

Table  I)  in  the  vicinity  of  0.8^7  Mev.   The  results  are  shown  in  Fig. 

11.  It  was  later  determined  that  a  much  thinner  neutron  target  was 

needed  to  resolve  the  structure  and  to  search  for  cusps  (see  Fig.  Ik; 

target  #5,  Table  I  was  used). 

Energy  and  Angle  Dependence  of  Threshold  Features 

We  discuss  here  the  experimental  results  as  they  pertain  to 
Wigner  and  "energy-averaged"  Wigner  cusps.  Frequent  reference  will  be 
made,  explicitly  and  implicitly,  to  the  predictions  of  cusp  theory, 
Chapter  II. 
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Fig.  10.  Neutron  differential  elastic  cross  sections  (CM)  of 
natural  lithium  versus  center  of  mass  (CM)  energy  for  five  angles  (CM) 
The  data  are  not  corrected  for  multiple  scattering.  The  original 
experimental  results  have  been  averaged  such  that  statistical  errors 
in  these  points  are  approximately  ±1$. 
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Fig.  11.  Neutron  differential  elastic  cross  sections  of  natural 
iron  versus  energy  for  four  angles.  The  data  is  not  corrected  for 
multiple  scattering.  Incident  neutron  energy  spread  is  approximately 
25  kev.  Statistical  errors  are  indicated  by  the  size  of  the  symbol. 
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Lithium  (92  .6#  II7) 

The  cusp  results  for  lithium  were  essentially  negative .  The 
theory  of  Appendix  B,  Case  I,  applies.  We  saw  in  the  preceding  section 
that  the  largest  percentage  cusp  effect  in  da  is  expected  at  forward 
angles.  Figure  12(b)  shows  that  an  experimental  effect  was  indeed 

present  at  3^°  (CM)  in  the  vicinity  of  0 A78  Mev,  the  first  excited 

7 
state  of  Li  •  However,  using  the  reaction  cross  section  magnitude  of 

29 
Freeman,  Lane,  and  Rose,   and  the  applicable  theory,  we  find  that  the 

observed  effect  was  at  least  a  factor  5  "too  large  to  be  explained  as 

a  Wigner  cusp;  an  effect  of  approximately  2  percent  is  predicted  50  kev 

from  threshold.  Figure  12(a)  compares  our  measurement  of  lithium  total 

27 
cross  section  with  that  from  Howerton's  compilation.    Although  the 

two  sets  of  data  agree,  no  Wigner  cusp  is  in  evidence  in  either.  We 
had  hoped  to  see  an  effect  in  the  total  cross  section,  because  the  re- 
sults of  Lane,  et  al.,  (see  Fig.  13)  appear  to  show  a  break  near  0„5 
Mev  in  the  coefficient  B  of  the  cross  section  expansion  (see  Chapter 
II).  The  energy  scale  change  in  Fig.  13  at  about  0.5  Mev  may  bias  this 
interpretation  of  their  data. 

Iron  (91^7$  Fe  ) 

With  iron  distinct  features  appear  in  all  cross  sections  at  the 

+  56 

energy  of  interest  (2  first  excited  state  Fe  ,  O08U7  Mev).  This  is 

exemplified  in  Fig.  Ik   where  a  resonance  feature  appears,  apparently 
at  an  identical  energy  (just  above  O.8U7  Mev)  in  the  total  cross  sec- 
tion, in  the  1^5°  differential  elastic  cross  section  and  in  the  reaction 
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Fig.  12.  Neutron  cross  sections  of  natural  lithium  versus 

energy,   (a)  Total  neutron  cross  section;  triangles  are  data  from 

27 

Howerton's  compilation   and  solid  circles  are  the  energy-averaged 

(A  E  =20  kev)  data  from  the  present  experiment,   (b)  3^°  CM  differ- 
ential elastic  cross  section  of  the  present  experiment. 
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Fig.  13.   Differential  cross  section  for  scattering  of  neutrons 

7 
by  Li  .   Reproduced  from  the  work  of  A.  Langsdorf,  Jr.,  R.  0.  Lane,  and 

20 

J .  E „  Monahan . 
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Fig.  Ik.      (a)  Neutron  differential  elastic  (1^5°)  cross  section 
of  natural  iron  versus  energy,  (b)  Total  neutron  cross  section  of 
natural  iron,   (c)  Reaction  cross  section  Fe   (n,  n1,  y  (.8*4-7  Mev))  Fe   , 
Statistical  errors  are  the  size  of  the  symbols  in  (a)  and  (b),  but  are 
±10$  in  (c).  The  smooth  curves  show  the  result  of  removing  a  calculated 
Wigner  cusp  effect  under  one  possible  assumption  on  the  phase  of  the 
scattering  amplitude. 
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cross  section  Fe  (n,  n')Fe   (0.847  Mev).  The  latter  cross  section 

10 
was  measured  by  Tucker,  et  al  ,  using  the  same  tritium  target.    An- 
other resonance  feature  which  seems  to  be  a  mirror  image  '  of  the 
aforementioned  resonance  across  threshold  is  evident  in  these  data, 
and  in  all  the  differential  cross  sections  at  other  angles  investi- 
gated [see  Fig.  15  (six  angles)] .  Application  of  the  most  general 
theory  is  hopeless  because  of  the  large  number  of  unknown  phases.  At 
least  six  phases  contribute  to  the  cusp,  even  with  a  restriction  on 
participating  incident  partial  waves  to  those  with  orbital  angular 
momentum  £     less  than  3;  see  Appendix  B»  A  reasonable  and  tractable 
special  case  of  the  theory,  which  might  apply  to  the  iron  data  is  dis- 
cussed and  explicitly  evaluated  in  Appendix  B,  Case  II.   Unfortunately, 
the  special  assumptions  involved  allow  cusp  analysis  only  over  a  nar- 
row energy  region,  entirely  above  threshold.  The  results  of  this 
analysis  are  superimposed  on  the  data  in  Fig.  15  for  differential  cross 
sections  (6  angles)  and  in  Fig.  l6  for  the  total  cross  section.   The 
results  are  neither  startling  nor  conclusive,  although  by  removal  of 
the  cusp  effect  the  peak  energies  and  widths  of  the  above  threshold 
resonance  in  the  differential  elastic,  total  and  inelastic  cross  sec- 
tions can  be  brought  into  exact  coincidence. 

Cerium-Barium- Zirconium 

We  consider  now  the  results  obtained  for  targets  of  higher 
atomic  number  in  the  region  where  the  optical  model  is  expected  to  ap- 
ply because  compound  resonances  are  closely  spaced- -cerium,  barium, 
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Fig.  15.   Neutron  differential  elastic  cross  sections  of 
natural  iron  versus  energy  for  six  angles.  The  data  is  not  corrected 
for  multiple  scattering,   (a)  gives  the  data  over  the  full  energy- 
range  of  the  run.  The  lined  box  shows  the  section  of  the  data  which 
is  expanded  into  (b).   In  (b)  statistical  errors  are  indicated  by 
the  symbol  size.   Open  triangles  show  the  result  of  removing  a  cal- 
culated Wigner  cusp  effect,  under  one  possible  assumption  on  the 
phase  of  the  scattering  amplitude. 
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Fig.  l6.   Total  cross  section  of  natural  iron.   The  size  of 
statistical  errors  is  indicated  by  the  size  of  the  symbol ■   Triangles 
show  the  result  of  removing  a  calculated  Wigner  cusp  effect  under  one 
possible  assumption  on  the  phase  of  the  scattering  amplitude. 
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and  zirconium.   The  theory  of  "energy- averaged  '  Wigner  cusps  is  then 
also  expected  to  apply  (see  Chapter  II). 

Composite  data  for  cerium  (88.5$  Ce   ),  given  in  Fig.  17,  is 
representative  of  the  cusp  analyses  we  present  for  all  three  elements . 
A  special  feature  of  this  figure,  however,  is  that  we  present  the 
original  data,  to  show  the  degree  of  scatter  on  interleaving  runs  (most 
of  which  is  not  statistical,  hut  experimentally  reproducible).   In  the 
remaining  figures  presented  for  this  group,  the  data  was  averaged  over 
an  energy  interval  suitable  to  reduce  statistical  errors  to  about  ±1 
percent  (refer  to  Chapter  IV ) .  This  procedure  standardized  data  presen- 
tation for  a  large  number  of  angles . 

This  and  other  figures  in  this  group  also  show  cross  section 

trends  predicted  by  the  optical  model.   These  were  obtained  from  the 

30 
work  of  Campbell,  et  al.    The  optical  model  shapes  were  interpolated 

from  very  few  points  given  in  the  reference  and  should  not  be  consid- 
ered a  precise  computation;  hence,  they  were  normalized  to  the  experi- 
mental data  at  threshold  energy.  The  magnitudes  of  the  predicted  cross 

30 
sections   were  consistent  with  the  estimated  experimental  errors 

(Chapter  IV )  and  the  corrections  given  in  Appendix  E„ 

A  cusp  effect  is  quite  apparent  in  Fig.  17(b)  by  comparison  of 
the  optical  model  trends  with  the  experimental  cross  section.  To  com- 
pare the  magnitude  of  the  cusp  effect  with  theory  one  could  subtract 
the  experimental  cross  section  from  the  optical  model  prediction  and 
plot  the  difference  as  a  function  of  neutron  energy.   But,  as  mentioned 
before,  the  optical  model  trends  are  somewhat  uncertain  and,  furthermore, 
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Fig .  17 •   The  experimental  inputs  to  a  cerium  "energy-averaged 
Wigner  cusp  analysis,   (a)  Energy- averaged  (A  E  =  ^40  kev)  experimental 
inelastic  cross  section  Ce   (n,  n1  y  (1.597  Mev))  Ce   .  '   (b)  Origi- 
nal data  for  elastic  95°  scattering  of  neutrons  from  natural  cerium. 
The  two  symbols  indicate  the  reproducibility  of  the  data  on  successive 
interleaving  energy  runs.  The  size  of  the  symbols  indicates  the  size 
of  the  statistical  errors,   (c)  Same  data  as  (b)  corrected,  point  by 
point,  for  a  calculated  'energy-averaged"  Wigner  cusp.   The  smooth 
curves  in  (b)  and  (c)  are  optical  model  predictions  of  Campbell,  et  al., 
normalized  at  1.6  Mev. 
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the  experimental  cross  section  has  natural  fluctuations,  which  will  be 
discussed  below.  Hence,  it  was  thought  more  appropriate  to  remove  from 
the  experimental  cross  sections  the  theoretical  cusp  effect  and  to  show 
(a)  that  in  the  resultant  curve  there  would  be  no  more  cusp  effect—at 
least  within  the  natural  fluctuations  of  the  cross  section,  and  (b) 
that  the  resultant  curve  would  follow  the  optical  model  trend,  which 
theoretically  represents  the  cross  section  without  threshold  effects. 

The  procedure  for  removing  a  theoretical  cusp  from  the  measured 
cross  section  was  to  multiply  the  reaction  cross  section  (i.e.,  inelas- 
tic neutron  scattering  cross  section)  by  the  isotopic  percentage  of  the 
isotope  of  interest,  by  the  appropriate  multiple  scattering  and  finite 
scatterer  correction  (the  reciprocal  of  Table  III  (Appendix  E)  values), 
and  by  the  theoretical  angle  dependent  factor  (Chapter  II) .   This  neg- 
ative quantity  was  subtracted  from  the  experimental  data  [see,  for 
example,  Eq.  (26)]. 

lUo 
Cerium  (88.5$  Ce   ) 

We  now  consider  the  specific  case  of  cerium,  for  which  Eq.  (26) 
applies.   Tucker,  et  al.,  measured  the  reaction  cross  section  [Fig. 
17(a)]  Ce   (n,  n')Ce    (l.6o  Mev),    As  previously  mentioned,  inci- 
dent d  waves  dominate  the  reaction  to  at  least  150  kev  above  threshold, 
hence  i '  =  2  in  Eq.  (26).  The  results  of  removing  the  theoretical 
cusp  effect  are  shown  in  Fig.  17(c).  There  is  little  doubt  that  the 
theoretical  "energy-averaged"  Wigner  cusp,  with  Zx    =   2  Is  a  satis- 
factory description  of  the  change  at  threshold.  Note  Eq.  (26)  predicts 
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that  if  incident  _s  waves  were  imagined  to  dominate  this  reaction,  rather 
than  the  actual  d  waves,  the  resulting  expected  cusp,  at  polar  angle 
90°,  would  he  a  factor  20/13  or  about  50  percent  larger.   The  plot  cor- 
responding to  Fig.  17(c),  corrected  for  an  expected  (i1  =  0)  cusp, 
would  show  an  upvard  break  at  threshold  and  would  follow  neither  the 
below  threshold  data  trend  nor  the  optical  model  trend.   Thus  we  can 
conclude  the  reaction  is  not  dominated  by  incident  _s  waves.   This  type 
of  information  can  be  useful  in  assigning  excited  state  spins,  for 
example . 

The  differential  elastic  cross  section  of  cerium  was  measured 
at  ten  angles,  versus  energy.   These  results  are  shown  in  Fig.  18. 
The  points  obtained  after  removing  a  calculated  cusp  are  also  given 
for  each  angle.   In  all  cases,  the  adjusted  data  so  obtained  follow 
the  trend  of  the  below  threshold  data  and  the  fluctuations  in  these 
is  the  same  above  and  below  threshold.   In  general,  the  optical  model 
curves  agree  nicely  with  the  modified  cross  sections;  that  is,  they 
follow  what  the  theory  defines  as  "cross  section  without  reaction." 
There  is  no  reason  to  expect  this  agreement  to  be  perfect,  since  the 
optical  model  was  constructed  to  explain  gross  structures  only,  over 
a  wide  range  of  atomic  number  and  energy.  We  should  recall  the  assump- 
tion, made  for  Eq.  (26),  that  the  fraction  (called  R  ,   in  Chapter  II ) 
of  the  reaction  probability  for  proceeding  from  compound  nucleus  states 
of  a  particular  J'   is  the  same  as  the  fraction  of  the  dominant  inci- 
dent partial  wave  states  (all  with  the  same  orbital  angular  momentum 
i ' )  which  have  J  =  J1.   Other  assumptions  for  RTi   lead  to  cusp 
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Fig.  18.   The  neutron  differential  elastic  cross  section  of 
natural  cerium  versus  energy  for  ten  angles.  The  data  is  not  corrected 
for  multiple  scattering.  Statistical  errors  are  ±1$.   The  points  are 
an  energy  average  (A  E  =50  kev)  of  the  original  experimental  data. 
For  each  angle  the  left  side  is  the  measured  cross  section.   The  right 
side  is  the  cross  section  with  a  calculated  "energy-averaged"  Wigner 

cusp  removed.   The  calculation  is  based  on  the  excitation  of  the  1.6 

11*0  +  ,        . 

Mev  state  in  Ce    .  No  cusp  contributions  from  the  0   (1-90  Mev;  or 

k     (2„08  Mev)  states  are  included.   The  smooth  curves  are  optical  . 

30 
model  predictions  of  Campbell,  et  al.,  normalized  at  1.6  Mev, 
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magnitude  predictions  differing  from  those  used  here  by  as  much  as  30 

percent,  at  certain  angles.  Although  not  the  case  in  the  present  work, 

it  is  possible  that  a  more  refined  experiment  could  be  used  to  measure 

R  ,;   such  a  measurement  could  assist  in  the  understanding  of  neutron 
J 

nucleus  interaction  mechanisms,  particularly  if  it  turned  out  that  a 
particular  reaction  preferred  one  value  of  compound  nucleus  J'  to 
another,  even  though  resonances  were  closely  spaced  and  possibly  over- 
lapping . 

Barium  (71.7$  Ba13  ) 

The  measured  cross  sections  of  barium  and  the  results  of  the 
appropriate  cusp  analysis  are  compiled  in  Fig.  19 .  The  reaction  cross 

1  qO  t  o Q-¥-  T  Q 

section  Ba   (n,  n')Ba    (1.427  Mev)  was  measured  by  Tucker,  et  al.  ' 

Again  the  spin  of  the  first  excited  1.^27  Mev  level  is  2  and  the  ground 

+ 
state  0  .  That  s   wave  neutrons  dominate  the  reaction  very  far  from 

threshold  is  not  so  clear  in  this  case,  since  the  reaction  shape  is 

distorted  by  resonance  structure.  However,  optical  model  calculations 

18 

show  that  barium  is  near  a  maximum  in  the  _s  wave  strength  function. 

As  in  cerium  then,  the  reaction  is  considered  to  proceed  by  the  &x    -   2 
incident  partial  wave.   The  results  of  theoretical  comparison  are  sat- 
isfactory. The  cusp  effect  is  smaller  than  in  cerium  because  of  the 
smaller  isotopic  percentage  and  the  smaller  reaction  cross  section  near 
threshold.  In  Fig.  19(c)  we  see  that  it  would  be  difficult  to  distin- 
guish a  cusp  at  all  were  it  not  for  a  comparison  with  the  optical  model 
trend  which  happens  to  be  rising. 


Fig.  19.  Neutron  differential  elastic  (90°  and  125°)  cross 
;ection£  of  natural  bari  nn  versus  energy,   (a)  Energy-averaged  (A  E  = 
20  :.(,  )  inelastic  cross  section  for  Ba   (n,  n'  r   (1.^27  Mev))  Ba  ^°. 
il  error  are  ±i°/o.     The  left  sides  of  (b)  and  (cv)  are  an 
a  era^.e  (A  E  -  2,  kev)  of  the  original  experimental  data.  This 
[   ot  corrected  for  multiple  scattering.  Statistical  erroi-o  are 
±1$.  The  right  side  of  (b)  and  (c)  show  the  results  of  removing 

calculated  "energy-averaged"  Wigner  cusp:;.  The  smooth  curves  arc  op- 

.       30 
al   -del  predictions  of  Campbell,  et  al.,  normalized  at  1.4-27  Mev. 
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Zirconium  (51.5#  Zr   i 

The  measured  cross  sections  of  zirconium  and  the  results  of 
cusp  analysis  are  compiled  in  Fig.  20.  The  inelastic  reaction 

Zr  (n,  n')Zr   (1-75  Mev)  was)  measured  by  Tucker,  et  al.  '   The  spin 

+  + 

of  the  1.75-Mev  first  excited  state  is  0  and  the  ground  state  0  . 

Tucker  also  shoved  that  p_  wave  outgoing  neutrons  dominate  the{  reaction 

not  far  above  threshold.  Hence  we  assume  that  the  reaction  proceeds 

90 
from  incident  p_  waves,  i.e.,   i'  =1.  For  Zr   the  low  isotopic  per- 
centage and  small  reaction  cross  section  near  threshold  combine  %o 
produce  a  small  calculated  cusp  effect.  At  60°  polar  angle  the  pre- 
dicted cusp  effect  is  too  small  to  be  distinguished  in  the  data.  Note 
also  that  at  60°  the  optical  model  trend  does  not  agree  well  with  the 
measured  cross  section  slope.  At  130°  the  cusp  is  distinguishable  and 
the  theoretical  cusp  fit  is  satisfactory.. 

Total  Cross  Sections 

The  theoretical  'energy-averaged"  Wigner  cusps  which  were  used 

to  explain  the  preceding  experimental  results  (Ce,  Ba,  Zr)  are  based 

22 

on  the  same  assumptions  that  underlie  the  optical  model*    These  as- 
sumptions also  lead  to  predictions  that  the  energy  average  of  the  total 
cross  sections  should  be  smooth  across  thresholds  [see  Eq.  (23)].   Our 
measurements  of  the  total  cross  sections  of  cerium  and  barium  are  given 
in  Fig.  21.  The  predicted  cusp  effect  for  the  total  cross  section  of 
zirconium- -even  if  the  energy- averaged  theory  would  not  apply--would 
be  so  small  that  we  did  not  investigate  this  case.  For  completeness, 
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Fig.  20.  Neutron  differential  elastic  (60°  and  130°)  cross 

sections  of  natural  zirconium  versus  energy,   (a)  is  a  1*0  kev  average 

90         90  19 

of  the  reaction  cross  section  Zr  (n,  n1)  Zr  (1.75  Mev) .  '   Statis- 
tical errors  are  ±5$-  The  left  sides  of  (b)  and  (c)  are  an  energy 
average  (A  E  =50  kev)  of  the  measured  cross  sections.  This  data  is 
not  corrected  for  multiple  scattering.  Statistical  errors  are  ±1$. 
The  right  sides  of  (b)  and  (c)  show  the  results  of  removing  calculated 

"energy-averaged"  Wigner  cusps.   The  smooth  curves  are  optical  model 

30 

predictions  of  Campbell,  et  al_,  normalized  at  1.75  Mev. 
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Fig.  21.   Neutron  total  cross  sections  for  natural  barium,  (a); 
cerium,  (b);  and  zirconium,  (c).  Statistical  errors  are  indicated  by 
the  size  of  the  symbols.   The  data  on  barium  and  cerium  were  measured 

in  the  present  experiment.   The  data  on  zirconium  is  taken  from 

27 

Howerton's  compilation.    In  (b)  circles  are  the  total  cross  section 

of  cerium;  triangles  are  integrated  elastic  cross  sections,  corrected 
for  multiple  scattering;  and  squares  are  the  sum  of  the  integrated 
elastic  and  of  the  reaction  cross  section  (corrected  for  isotopic  per- 
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27 
the  zirconium  total  cross  section  data  from  Howerton   is  also  presented 

in  Fig.  21,  showing  the  inadequacy  of  widely  spaced  energy  points  for 
the  present  study.  The  data  for  barium  and  cerium  were  taken  with  an 
energy  spread  of  about  30  kev,  yet  they  are  certainly  not  smooth  func- 
tions of  energy,  displaying  marked  structure .  The  mean  values  of  these 
fluctuations  do  not  seem  to  drop  markedly  as  threshold  is  crossed,  sug- 
gesting agreement  with  the  theory ■   If,  however,  the  theory  of  energy- 
averaged  cusps  does  not  apply,  but  one  assumes  instead  a  pure  Wigner 
cusp  phenomenon,  the  cusp  predicted  is  at  most  a  5  percent  effect  which 
would  be  masked  by  the  fluctuations.  Hence,  the  total  cross  section 
data  cannot  be  considered  a  verification  of  either  assumption,, 

That  the  fluctuations  are  not  experimental  is  supported  by  the 
differential  cross  section  data  for  cerium  at  forward  angles  (Fig.  18) 
which  shows  the  same  structure.  Note  that  the  structure  does  not  con- 
tain recognizable  resonance  shapes.  It  is  possible  that  the  consider- 
ations of  Eric son  on  statistical  fluctuations  in  cross  sections  might 

31 
apply  here .    He  shows  that  if  a  compound  system  is  characterized  by 

many  overlapping  resonances,  each  of  width  r,   and  with  relative  ran- 

21 
dom  phase  of  the  value  quantities  y   of  resonance  theory,   cross 

A, 

sections  would  have  a  fluctuating  structure  of  average  width,   P. 

In  the  case  of  cerium  only,  sufficient  differential  elastic 
data  was  taken  to  permit  integration  over  all  solid  angle.   The  total 
elastic  cross  sections  so  obtained  for  several  energies,  after  correc- 
tion for  multiple  scattering,  are  plotted  in  Fig.  21(b).  The  drop  at 
threshold  (1.6  Mev)  is  evident,  as,  in  accord  with  any  theory  of 
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threshold  effects,  flux  conservation  alone  requires  the  total  elastic 

cross  section  to  drop  above  threshold  for  a  reaction.   In  the  case  of 

1^0 
Ce   ,  the  experimental  differential  cross  sections  all  drop,  so  the 

total  elastic  cross  section  does  also.   It  is  reassuring  that,  when 
the  inelastic  scattering  cross  section  is  corrected  for  isotopic  per- 
centage and  added  to  the  integrated  elastic  cross  section,  the  result- 
ant cross  section  (i.e.,  the  total  cross  section)  is  smooth  across 
threshold.  The  agreement  in  magnitude  between  the  latter  curve  and 
the  directly  measured  total  cross  section  is  considered  somewhat 
fortuitous . 


CHAPTER  VI 
CONCLUSIONS 

The  experiments  have  shown  that  under  suitable  conditions  the 
threshold  effects  expected  in  neutron  elastic  scattering  can  be  ob- 
served. From  an  experimental  point  of  view  particular  attention  must 
be  paid  to  the  choice  of  energy  resolution  consistent  with  observing 
a  cusp  effect  in  reasonable  time. 

Explicit  theoretical  expressions  have  been  obtained  for  ex- 
pected Wigner  and  "energy-averaged"  Wigner  cusp  effects  under  general 
angular  momentum  conditions,  within  some  reasonable  restrictions  imposed 
in  the  derivation.   The  effects  were  found  to  be  contained  in  the  coef- 
ficients,  BT,   of  the  Blatt-Biedenharn  formulation  of  the  general 
scattering  problem.   Since  many  survey  experiments  involving  neutron 
scattering  have  been  analyzed  in  terms  of  the  BT   coefficients,  we 
make  the  obvious  suggestion  that  such  results  be  examined  near  excited 
state  energies  for  possible  tests  of  cusp  theory. 

The  expected  Wigner  cusp  in  lithium  cross  sections  was  expected 
to  be  an  approximately  2  percent  effect,  50  kev  from  threshold.   Our 
experiment  was  not  sufficiently  refined  to  pick  out  an  effect  of  this 
size. 

Resonance  structure  in  the  measured  cross  sections  of  iron  made 
a  Wigner  cusp  interpretation  ambiguous.  Features  occurred  at  the  thresh- 
old energy,  but  their  analysis  in  terms  of  cusp  theory  was  inconclusive. 
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The  measured  differential  elastic  cross  sections  of  Zr,  Ba, 
and  Ce  dropped  at  inelastic  thresholds.  Analysis  of  the  drops  in 
terms  of  "energy-averaged"  Wigner  cusp  theory  resulted  in  good  agree- 
ment between  theory  and  experiment „  This  indicates  that  the  usual 

energy -averaging  assumptions  underlying  the  optical  model  are  indeed 

7 

applicable  to  these  nuclei.  Also,  the  cusp  theory  shows  that  there 

cannot  be  a  large  amount  of  direct  interaction  between  incident  neu- 
trons and  these  nuclei  at  the  neutron  energies  used  in  these  experi- 
ments (^  2  Mev),  or  else  below  threshold  effects  would  also  be  seen. 


APPENDIX  A 

THE  PHYSICAL  BASIS  FOR  THE  S-MATRIX  AND  ITS  USE 
IN  CROSS  SECTION  FORMULAS 

We  will  review  the  principles  of  scattering  analysis  and  the 

12 
definition  of  the  S-matrix  as  presented  by  Blatt  and  Biedenharn. 

They  pointed  out  that  general  expressions  for  cross  sections  in  terms 

of  the  S-matrix  are  rarely  given  in  standard  texts . 

The  general  formulas  are  applicable  to  any  collision  process 
in  which  two  particles  collide  and  two  particles  emerge.  However,  in 
this  appendix,  and  after  a  certain  stage  in  the  development,  only  the 
special  case  of  elastic  scattering  with  no  i  change  will  be  con- 
sidered. 

The  channel  index  a.     defines  the  type  of  incoming  particle  and 
the  state  of  the  struck  nucleus.   The  channel  spin,   s,   is  the  total 
spin  angular  momentum  of  the  channel,  obtained  by  vector  addition  of 
the  intrinsic  spin,   i,   of  the  incident  particle  and  the  spin,   I, 
of  the  struck  nucleus.   The  orbital  angular  momentum  (in  the  center  of 
mass  system),   i,   is  combined  with  s  to  form  the  total  angular 
momentum,   J.   Only  in  this  particular  appendix,  primes  will  be  used 
to  refer  to  outgoing  or  final  state  quantities.   J  is  conserved  during 
the  collision,  hence  the  values  of  i '   and  s'  must  be  consistent 
with  J'  =  J.   Similarly,  the  parity  jt,   of  the  initial  system,  is 
conserved,  and  n    =  it '  . 
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The  probability  amplitude  S  ,    ,  0  t    0     for  a  collision  with 
total  angular  momentum  J,   from  channel  OLsZ     to  channel  a's'i'   is 
referred  to  as  an  element  of  the  scattering  matrix.   The  scattering 
matrix,  or  S-matrix,  is  defined  by  considering  asymptotic  forms  of  the 
collision  wave  functions,  as  follows. 

Consider  one  definite  value  of  J  and  J  =  M.  At  any  particu- 
lar energy,   E,   the  channel  wave  number,  k  ,   and  relative  speed, 
v  ,   in  each  channel  a  are  given  by  the  energetics  of  the  reaction. 
Only  those  channels  which  are  "open"  (k  real)  are  considered,   Let 
cp    be  the  product  wave  function  of  the  target  nucleus  and  incident 
particle.   The  spin  and  angle  dependence  of  a  wave  function  of  total 
angular  momentum  J,   J  =  M,   orbital  angular  momentum  Z,      and  chan- 

Zl 

nel  spin  s  is  given  by 


*"* = .  ii .  E-s {'sm*  m^s  m)  We'v)  x*'ms  (A1) 

Z  s 

where  X      is  the  incident  spin  wave  function  and  Y.   is  the  well 
s,m  *  £m 

s 

known  spherical  harmonic  of  angular  momentum  Z     and  projection  m. 
Using  these  definitions,  the  most  general  wave  function  in  channel  OL,  s 
with  total  angular  momentum  quantum  numbers  J,  M,   consists  of  the 
superposition  of  incoming  and  outgoing  spherical  waves .  Asymptotically 
(at  sufficiently  large  distances)  we  can  write  (for  neutral  particles) 


*   (jM)  = L-r   <aM   ^ 

aS       r  (v  )*  "  J  s  as 


a  a 


x  iAas  lM  ^l-'Va  "  5  '  «»  -  Bas  ™  exp[+i(Va  - 

\  I  «)]}     (A2) 


Ill 


Substitution  of  Eq.  (A2)  into  the  wave  equation  yields  a  unique  connec- 
tion between  the  outgoing  wave  coefficients,  B,  with  the  ingoing  wave 
coefficients,   A.   Thus 

BaV?  =lZL,    saVr;cJAasf  <A3) 

a  s  £  ' 

Eq.  (A3)  defines  the  scattering  matrix,   S.  With  these  definitions  it 
is  evident  that  S  is  unitary  (conserving  the  wave  function  normaliza- 
tion in  time),  symmetric  (reciprocity  theorem),  and  independent  of  M. 
It  remains  to  decompose  the  incident  plane  wave  into  the  func- 
tional forms  we  are  using.   The  incident  wave  is  cp   exp(ikz)X 

lxo  s  *  m 

s 

First,    one  writes 

s         i!=0  *       7      '  s 

where 

i+s      J 


Y,  ft(6)X     =  Z       Z        (isom  lis  JM)UT, 
i,0    s,m„    T_, „   „,  „_  T      s1       Jjis 


Using  the  asymptotic  form  of  J   i   for  kr  >  >  £     the  following  ex- 
pression  is  obtained  for  the  plane  wave  in  channel  a, s  with  spin 


direction  given  by  m  ; 


1 


oo  J  J+S 


exp(ik   z    )X         cp       s    iL*L  cp 
■^v     a  a    s,m  ras        k  r      Yas   _  _  .,     _    .   ,  _ 

'    s  a  a  J=0  M=-J  £=  J-s 


X  (isom   |is  JM)   i    (2i  +  l)2  U       M  (AU) 

S  0  J  iJs 

x|exp[-i(kara  -iiy)]-   exp[+i(karQ;  -  \  i  «)]| 


112 


Comparing  Eq.  (A^)  with  the  standard  form,  Eq.  (A2),  we  see  that  the 
amplitudes  of  the  ingoing  spherical  waves  are  given  by 

A  f  =  i)c  (it  v  )2  (isom  lis  JM)  ii(2i  +  l)2  (A5) 

asi     a   a      s1 

The  amplitudes  of  the  outgoing  spherical  waves  are  determined  from  Eq. 
(A3)-  As  in  the  main  text,  we  now  specialize  to  elastic  scattering 
with  no  change  of  £     or  s .   They  are 

™  X   J+s 

B  f   =  i  *  (*  v  )2  2 

asi     ocK      or  .  i  T  i 

i=| J-s| 

fl         4  T 

X  (isom  is  JM)  i  (2i  +  l)2  S  „        „  (A6) 

v    s1     y   v     '   asi;asi  v  ' 

In  order  to  investigate  cross  sections  we  only  need  that  part  of  the 
plane  wave  solution,  Eq.  (kh) ,   which  is  due  to  the  scattering  interac- 
tion: i.e.,  we  want  to  obtain  \If        ,  from  if  =  it .      .,  + 
'                       T scattered           incident 

•it  ..   Therefore,  we  substitute  Eqs .  (A5)  and  (A6)  into  Eq.  (A2) 

scattered 

and  subtract  the  incident  wave,  Eq.  (kk) ;   the  result  is 


1  oo    j    J+s 

2  J- 
..  CD 

scattered      a  r  ^as 


*...*+—*  =  i  Kv  *2  ~  <P~.  s 


"a    j=o  m=-j  i=|j-s| 


X  (isom  is  JM)  i  (2i  +  l)2 

s ' 


X  e«p[i(V«  -  |  *  «)1  X  [1  -  S^    ]^         (A7) 


M 


If  the  spin  angle  function  l^-r^   is  decomposed  using  Eq.  (Al),  then 
Eq.  (A7)  can  be  used  to  define  the  scattered  amplitude,   a    r. 


Ctsm 
s'   s 


*   (asm1)  =  I*  exPlkr    cp 
Y  sc    s      a   r      as 


X   I   a   ,     (e  q>)  X   ,  (a8) 

TKsm  :asm  v  /   s ,  m 
m  =-s     s'    s        '  s 
s 
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where  q  is  explicitly 

00  J   J+s      i 

a    ,     (0,cp)  =  E   E    E       E 
Xism-  :asm   ,Y/    T_r.   M   T  -  iT   i     , 
5'    s        J=0  M=-J  x=|  J-s  |  |J.=-i 

1  l 

X  7^(2^  +  l)2(isom  |  is  JM)(isMm*Jis  JM) 

S  o 

*  (1  -  ScJjasi)  V'*5  U9) 

The  differential  cross  section  for  the  scattering  asm  ->asm'   is 

s      s 


given  by 

da   i      =  k 
asm  :asm     a 
s'   s 


a    ,     (6,q>)   &Q  (A10) 

Xtsm  ;0"" 


:asm 
s'   s 


For  unpolarized  incident  beams,  we  must  average  Eq„  (AIO)  over  initial 

spin  direction  m   and  sum  over  final  spin  direction  mT: 
s  s 

,    s     s 

do      =  (2s  +  1)     E     Z       do        ,  n  (All) 

as;Os  t      asm  : asm 

*  m  =-s  m  =-s      s'   s 

s     s 

This  expression  is  greatly  simplified  by  sum  rules  for  Clebsch-Gordan 
coefficients,  and  it  is  from  this  calculation,  as  performed  by  Blatt 
and  Biedenharn,  that  one  obtains  the  useful  final  cross  section  expres- 
sions, involving  Z  coefficients. 

Finally,  to  obtain  the  differential  elastic  cross  section,  one 
must  average  over  the  possible  values  of  incident  channel  spin,   s. 

1+1        2s  +  1 

daa;a  =  Z|T  .,  (21  +  l)(2i  +  l)  daas;as         (A12) 
'    s  =  |I-i|       'v     '     ' 

Equations  (A12)  and  (All),  as  explicitly  worked  out  in  Ref .  12,  are 
just  Eqs.  (l),  (2)  and  (3)  of  the  main  text. 


APPENDIX  B 
SPECIAL  CASES  OF  INTEREST  TO  THE  PRESENT  EXPERIMENT 

Case  I.  Cusp  effects  in  neutron  scattering  from  Li 

near  the  threshold  for  excitation  of  the  first 

n 
excited  (O.U78  Mev)  l/2  state  of  Li  (ground 

state  spin  and  parity,  3/2  ) . 

The  large  resonance  in  the  elastic  cross  section  at  about  0.25- 

+ 
Mev  neutron  energy  has  been  analyzed  and  found  to  correspond  to  a  3 

state  in  Li  .    The  reaction  cross  section  Li  (n,  n')Li   (O.V78  Mev) 

has  a  broad  resonance  at  1.35-Mev  neutron  energy,  probably  associated 

+       2Q 
with  a  compound  1  state.    Also,  the  inelastic  scattering  is  domi- 
nated near  threshold  by  outgoing  s   waves,  with  compound  angular  momen- 

JL  -  2Q 

turn  J  =  1  .    Finally,  any  process  involving  incident  partial  waves 

with  £  >   1  is  unlikely  at  0.5  Mev  because  kR  «  0  3-  We  conclude 
that  in  the  vicinity  of  0A8-Mev  incident  (center  of  mass)  neutron 
energy  only  s   and  p  waves  participate  significantly  in  elastic  scatter- 
ing. Also,  incident  _s  waves  dominate  the  inelastic  scattering.   The 

latter  conclusion  results  from  the  known  spins  and  parities  of  the 

7 
ground  and  first  excited  state  of  Li  ,  and  the  assumption  that  no 

process  with  i  >  1  is  likely .  With  _s  waves  dominant  in  the  outgoing 
reaction  channel  the  only  possible  incident  channel  leading  to  the  re- 
action is  characterized  by  J'  =  1,  s'  =  1,  and  £'  =  0.  Using  these 
values,  and  I  =  3/2,  i  =  l/2,   in  Eq.  (12),  one  obtains 


111+ 
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A  B  (1;1)  =  -  ^  Z  Z2(i  J  0  1,  1  L) 


.1    ,J 


sin(2501  -  5^  . 
*  2  sin  5L  *r01  {    f  1    5J  }  (B1) 

Restricting  ourselves  to  i  <  2,   as  indicated  above,  we  need  these 
Z  coefficients 

Z2(0101,10)  =  3,      Z2(  1101,11)  =  3,   Z2(  1201,11)  =  5 

and  all  other  Z  (iJ01,lL)  values  are  zero.   Therefore, 

2  sin  om 


'V1;1)  =-I4i{.  9i    } 


sin  250, 


iin(2501  -  b1±\ 


i    r        U1     L1-  i 


01   ir 


sin(26^  -  62  ). 

5      2   f      01  11 

+  3  Sin6nL  ,p*1  ^ 

W(2S   -  6   ) 


(B2) 


Equations  (B2)  might  be  compared  to  data  such  as  that  of  Langsdorf, 

20 
Lane  and  Monahan,   reproduced  herein  as  Fig.  13.   Unfortunately,  the 

energy  scale  in  this  figure  changes  in  the  region  of  interest  (0.5  Mev 

lab).  We  present  the  figure  primarily  to  illustrate  the  utility  of 

cusp  expressions  in  the  form  of  Eqs.  (B2).   The  cusp  prediction  for 

the  differential  elastic  cross  section  is  obtained  by  substituting  Eq. 

(B2)  in  Eqs.  (9-H)  to  obtain 
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Ada--(3/8)((arW))  expM«)      {     .""I10} 


2   sin     6. 


+  2  P  (cos  0) 


(5/3)   sin  6 


10 


sin(25oi  "  6ir 
11 


1  K2 

oi  -  811- 


+  sin  5 


sin(25io  -  sii} 


cos(2510  -  &1;L; 


cos(25~. 


dfl 


(B3) 


Using  Eq.  (l6),  we  obtain  the  total  cross  section  cusp  term, 

.1 

'01 


cos  26' 


f  cos  ^m       ~] 
AaT  =  (3/8)(0r^)exp{      °\  } 

*  v  -  sin  26^-,  J 


where   (a  „  ^)     in  Eqs .  (B3)  and  (B^)  is  the  experimental  reaction 
ri=0  exp 

cross  section  near  threshold. 

7 
For  the  case  of  Li  we  recall  the  restriction  imposed  on  the 

derivation  leading  to  Eq.  (12),  that  the  amount  of  scattering  with 
change  of  s,  hut  no  change  in  energy  or  orbital  angular  momentum  &, 
be  small.   There  are  no  good  reasons  for  believing  this  to  be  true  in 
the  present  case.  Change  of  s  =  I  +  i  from   |s|  =  2  to   |s  |  =  1 
(s  means  outgoing  "elastic"  channel  spin),  or  vice  versa,  in  the  Li 
case  is  essentially  a  spin  flip  followed  by  a  recoupling  of  s  and 
Z{  £   ^  l),   and  this  is  just  the  type  of  interaction  involved  in  the 
reaction  leading  to  the  first  excited  state.  We  have  seen,  however, 
that  in  fact  s   waves  (i1  =  0)  dominate  the  Li  reaction;  thus  "spin  flip" 
elastic  scattering  (which  requires  £   ^  l)  does  not  couple  to  the  reac- 
tion. Hence,  Eqs.  (B3)  and  (Bk)   apply.  Using  the  reaction  cross  section 


117 


29 
value  (100  kev  above  threshold)  of  Freeman,  lane,  and  Rose,   and  the 

present  experimental  values  of  da  (3^°  CM),   we  can  estimate  the  ex- 
pected percent  effect  to  be 

A  da/da°  =  (3/32  jt)  a  /da°  =   (3/32  it)  .05/ .05  =   3  percent  effect 

(~  100  kev  from  threshold) 
Correspondingly,  for  the  expected  cusp  in  a   we  obtain 


A  aT/aT°  -  (3/8)( -05/1-1)  *  2  percent  effect 

(100  kev  from  threshold) 


Case  II .  Cusp  effects  in  neutron  scattering  from 
iron  near  the  threshold  for  excitation  of  the 
first  excited  state  (2+,  O.8V7  Mev)  of  Fe5 
(ground  state  spin  and  parity  0  ) . 

Equations  (9)  to  (12)  for  pure  Wigner  cusp  effects  apply  to 
this  case  because  the  resonance  structure  in  the  cross  sections  is  ap- 
parently resolved.  We  must  compare  Eq„  (9);  the  expression  for  change 
in  cross  section,  with  experiment.  There  are  too  many  possible  phases 
involved  in  a  general  treatment  of  Eq.  (9)  to  make  cusp  interpretation 
meaningful;  therefore,  we  sought  a  tractable  and  reasonable  simplifica- 
tion. 

The  experimental  data,  see  Fig.  Ik,    display  a  resonance  at  an 
identical  energy  in  the  reaction  cross  section,  in  the  total  cross  sec- 
tion, and  at  six  angles  in  the  differential  elastic  cross  section.   On 
this  evidence,  we  assumed  that  in  the  vicinity  of  this  resonance  the 
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scattering  and  reaction  proceed  through  the  same  compound  state  of  total 
momentum,   J1,   and  are  dominated  by  one  incident  partial  wave,  Z\ 
With  this  assumption,  and  substituting  I  =  0,  i  =  1/2,  s'  =  1/2  in 

Eq.  (12)  we  obtain 

o      .  2  ft  J' 

2  sin  S£,i         y1 

A  B  (a  i;  a  i)  =  -  t  Z2(i'  J'  i'  J';  J  L)  4        Jt  [o"<l'A  (B5) 

L  sin  2  6  .  ,i_  J 

^  2 

where  J'  must  be  either  i'  +  ^  or  i '  -  \,      depending  on  the  nature 
of  the  compound  resonance.  By  substituting  Eq.  (B5)  in  Eqs .  (9)  to  (ll), 
one  obtains 


A  da 


2   J' 
/  J'   n      2  sin  5.ii 

(arTiW   r        J* 


a;a  "    kit  \    .  ORJ' 

'  ^   sin  2b e ti 

*  2 


2  Z  U'  fgj^'/i)  *  L)  PL(cOS  ^dQ  (B6) 


X 
L=0 


j" 
where  (cr  „,i)     is  the  measured  reaction  cross  section, 
ri  2  exP 

The  optical  model  indicates  a  maximum  in  the  s  wave  strength 

r   18  19 

function  at  A  w  5°>    ana-  "the  analysis  of  Tucker,  et  al.,    shows 

that  s_  wave  dominance  of  the  outgoing  neutrons  is  consistent  with  the 

measured  inelastic  cross  section  in  Fe   .  Knowing  that  the  excited 

state  spin  is  2  ,  we  conclude  &x    =   2  and  J1   is  either  5/2  or 

3/2.  For  J'  =  3/2; 

( o   3/2  j     2  sin  s|i2 

A  \a  =  "       ^f  ^  {  s.n  253(2  }  ^  +  p2  <cos  0)]da  (S^) 

2"2 
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For     J'    =  5/2: 


A  da. 


O         C  If 

2   sin     6pr 
Url'=2jexp   r 


a;  a 


5/2 


Jin  26^ 


2f 


1+|  P2(cos   0)    +  |  Pu(cos   9) 

(B7b) 


Using  Eq.  (l6),  we  obtain  the  expected  cusp  effects  in  the  total  cross 
section.  For  J1  =3/2; 


AaT  ^riUW 


cos  25 


-  si 


3/2 
in  26^ 


(B8a) 


For  J 


5/2: 


A  am  = 


cos  26^{2 

5/2       f      ^ 

^ri'^expl    .   OR5/2 
^  V.  -  sin  25^1 

22 


(B8b) 


Since  we  have  assumed  isolated  resonance  behavior  in  the  cross  sections 
in  this  case,  we  can  write  the  phases  5  as  explicit  functions  of  energy, 
Defining  e.,-1  =  (E  -  E  )/(rx,/2),  where  E   is  the  resonance  energy 

0  j    a   oJ  o 

and  rTi   is  the  full  width  at  half  maximum,  we  find 

2  sin2  5J'  =  2/(1+  e2)  , 

cos  2  6J;'   =  -  (1  -  €2)/(l  +  e2), 

T*  2 

-  sin  2  6   =  2  e/(l  +  e  )  , 


(B9) 


where  e  ranges  from  -  »  to  +  °°  across  the  resonance. 

Below-threshold  effects  are  possible  in  both  the  differential 
elastic  and  total  cross  sections  if  the  special  assumptions  are  valid. 
Unfortunately,  this  cannot  be  checked  in  the  present  experiment,  since 
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another  resonance  feature  appears  just  below  threshold,  whose  spin  is 
unknown . 

The  above-threshold  feature  can  be  analyzed  for  possible  cusp 
effects  within  the  framework  of  these  assumptions,  by  using  Eqs .  (B7) 
with  Eqs.  (B9)  substituted. 

A  da 


1  +  e  (BIO) 


X 


r(l  +  P2(cos  0),     if     J*    =  3/2 


L(i  +  (8/7)  P2(cos  6)   +  (6/7)  P^cos  9)),    if  J'    =  5/2  J 


dft 


^°T  =  (°ri^U^77f  (B11) 

Equations  (BIO)  and  (Bll)  show  an  interesting  property.   The 
cusp  term  in  expressions  for  elastic  cross  sections  (differential  and 
integrated)  is  negative  over  the  entire  resonance.   The  total  cross 
section  expression  is  negative  for  -  1  <  e  <  1,   and  positive  for 
e  <  -  1  and  -  e  >  1.   This  can  be  interpreted  simply  as  an  increase 
in  the  total  cross  section  resonance  width  due  to  the  reaction.   How- 
ever, the  experimental  neutron  energy  spread  («  7  kev)  was  too  broad 
to  allow  measurement  of  such  a  difference  in  the  widths  of  the  elastic 
and  total  cross  section  resonance  features. 

In  conclusion,  we  remark  that  the  below-threshold  feature  (see 
Fig.  Ik,    for  example)  which  we  have  described  as  "mirrored"  across 
threshold  is  just  the  kind  of  cusp  effect  one  might  expect  from  the 
theory.  We  were  disappointed  that  our  cusp  analysis  did  not  demonstrate 
this  feature  to  be  a  cusp  anomaly. 


APPENDIX  C 


A  PROOF 


We  wish  to  show  that  for  the  case  of  many  resonances  in  the 
averaging  interval  one  finds  under  certain  conditions 
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where  6   is  the  Kronecker  delta.  Expanding  the  left  side  of  (Cl), 
we  obtain 


j»       j  r  sin(28^,  -  6  ) 
(a      ,  2  sin  8  \  , 

ri        J   Lcos(26^  -  8^) 


=  <ari' 


J '       J         J '      2  J 

sin  26-i  sin  52  -  cos  28  ,  2  sin  6. 

J '       J         J '      2  J 

cos  25.1  sin  26  +  sin  28.,  2  sin  8. 


(C2) 


For  simplicity,  let 


Asa 


sin  28    , 
J'  J  l 


B  ■  a 


J' 
ri* 


-   cos  28., 
sin  28  . , 


a  a  sin  28 


J 


d  =  2   sm     8. 


121 


122 


Then  expression  (C2)  becomes 

=  (  Aa  +  Bb  >  (C3) 

We  now  assume,  for  £  £  £'     and/or  J  /  J1   and  within  an  averaging 
interval,  that  the  quantities  A,  B,  a,   and  b  may  be  written  as  the 
sum  of  a  mean  value  and  a  fluctuating  part;  for  example,   A  =  A°  +  A  A, 
where   (  A  A  )  =0.  Thus,  expression  (C3)  may  be  written 

(   A°a°   +  B°b°   +  A°A.a  +  a°A  A  +  B°A  b  +  b°A  B  +  AAAa+ABAb)    {ck) 

Meyerhof  showed  in  Ref .  7  that 

,  f  -  cos26  n 
(oJn,\  TIJ  [)  m ■<  B  >«0  (20) 


Also 


^   cos  25 . t 

on  the  same  basis  as  the  proof  of  (20).  Therefore,   A0  =  B°  «  0  and 
expression  (CU)  is  reduced  to 

(AAAa  +  ABAb)  (C5) 

where  we  have  also  used  (  a°A  A  )  =  (  b°A  B  )  =  0. 

A  A  and  A  B  are  functions  of  5  ,  while  A  a  and  A  b  refer  only 
to  6  .  For  £   /  £x      and/or  J  jt  J1,   the  fluctuation  over  the  averag- 
ing interval  in  functions  of  5.  are  uncorrelated  to  those  in  func- 

JO 
J'  /    N 

tions  of  6.,;  thus,  expression  (C5J  is  equal  to  zero,  and  we  have 


sin  28  , 


,  (■   sin  25  , 


proved  Eq.  (Cl)  for  this  case 
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For  J  =  J'   and  £   =  £'      the  LHS  of  Eq.  (Cl)  reduces  to 

.  2  .J* 
r£ 


<<^:{2sln  v}>, 


sin  2  bJn, 


and  ve  may  use  the  result  proved  in  Ref .  7> 

2  J1 
2  sin  6.,         /   J'\ 

r*   L  sin  2  6 , ,  J      L    0 


which  completes  the  proof  of  Eq.  (Cl). 


APPENDIX  D 

GAMMA  RAY  SUPPRESSION  PREAMPLIFIER 

A  space  charge  saturation  gamma-ray  suppression  preamplifier 
was  used  to  distinguish  neutron  induced  from  gamma-ray  induced  pulses 

in  the  detector.   The  circuit  is  shown  in  Fig.  22.   It  is  basically 

32 
that  of  Owens,   modified  to  include  more  external  controls,  a  diode 

clipper  on  the  nonlinear  output,  and  an  integrating  network  for  the 

remaining  positive  overshoot.   The  external  control  on  resistor  R 

33 
was  suggested  by  the  work  of  Funsten  and  Cobb. 

Certain  organic  scintillators  emit  light  under  particle  bom- 
bardment, which  is  composed  of  an  initial  spike  of  fluorescence  with 

about  50m|a  sec  decay  time,  followed  by  one  or  more  slower  components 

3k 
of  greater  than  200  mu-  sec  decay  times.    Ion  recombination  may  ac- 

35 
count  for  the  slow  components .    The  ratio  of  the  amplitude  of  the 

initial  spike  to  that  of  the  first  slow  component  is  a  function  of  the 
type  of  particle  causing  the  fluorescence . 

The  circuit  in  Fig.  22  is  designed  to  produce  linear  voltage 
pulses,  i.e.,  voltage  pulses  which  are  proportional  to  the  amplitude 
of  the  initial  fluorescence,  at  dynode  10.   The  remainder  of  the  re- 
sistor chain  is  arranged  deliberately  to  produce  space-charge  satura- 
tion in  the  region  between  dynode  Ik   and  the  anode,  during  the  initial 
spike.  When  the  slow-component  cascade  electrons  arrive  in  the  region 
of  space-charge  saturation  they  are  sufficiently  repelled  to  produce  a 

12k 
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Fig.  22.  Circuit  diagram  of  a  space  charge  saturation  gamma 
ray  suppression  preamplifier.   The  four  externally  controlled  elements 
are  marked  "FOC,"  "ACC,"  "R  , "  and  "DYN  Ik."      "LIN"  means  linear 
pulse  output.    y  supp"  means  the  gamma  ray  suppression  (nonlinear) 
output . 
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positive  overshoot  pulse  at  dynode  Ik.      Thus,  the  height  of  the  positive 
overshoot  is  a  measure  of  the  ratio  of  slow  component  to  initial  spike 
amplitude  and  it  is  possible  to  use  this  effect  to  distinguish  the  type 
of  particle  which  initiated  the  fluorescence. 

One  is  particularly  interested  in  distinguishing  neutron- 
initiated  pulses  (i.e.,  from  recoil  protons)  from  pulses  caused  by 
gamma-rays  (i.e.,  due  to  electrons),  and  the  method  has  proved  useful 
for  this  purpose.   Our  particular  circuit  is  arranged  such  that  the 
negative  pulse  at  dynode  Ik,    caused  by  the  initial  spike  for  either  in- 
cident neutrons  or  gamma- rays,  is  clipped  by  a  diode.  The  remaining 
positive  overshoot  is,  of  course,  no  longer  a  linear  effect,  and  we 
have  chosen  to  integrate  this  pulse  in  an  attempt  to  maximize  the  dis- 
tinction between  neutron  and  gamma-ray  initiated  pulses . 

The  circuit  was  effective  in  gamma  suppression  for  neutron 
energies  greater  than  about  0.3  Mev;  however,  some  care  must  be  taken 
in  biassing  out  gamma-ray  induced  pulses  and  to  select  neutron  caused 
pulses,  particularly  when  working  at  low  neutron  energies. 

An  example  of  the  nonlinear  spectrum  obtained  at  1.7-Mev  neu- 
tron energy  is  given  in  Fig.  23.  A  nonlinear  spectrum  of  gamma- ray 
initiated  pulses  is  also  shown  for  comparison.   These  gamma-rays  were 
produced  under  Van  de  Graaff  machine  operating  conditions  with  the 
proton  beam  incident  upon  a  stopper  in  front  of  the  tritium  target. 

While  working  at  low  neutron  energies  (<0.5  Mev)  we  found  that 
very  small  dynode  Ik   voltages  were  necessary  to  provide  good  gamma-ray 
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Fig.  23.   Nonlinear  pulse  spectra  obtained  at  E  =1.7  Mev 
from  a  space  charge  saturation  gamma  ray  suppression  preamplifier. 
The  detector  crystal  was  in  the  0°  neutron  beam  direction.   The 
solid  circles  show  the  nonlinear  spectrum  with  the  proton  beam  on 
the  T  target  (T(p,n)He  ).   Open  circles  show  the  nonlinear  spectrum 

with  the  proton  beam  blocked  by  a  tantalum  beam  stopper  in  front  of 

3 
the  T  target. 
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suppression.  This  produced  smaller  nonlinear  output  pulses,  which 
could  be  compensated  by  increasing  the  high  voltage  on  the  photomulti- 
plier,  but  which  eventually  represented  an  essential  limitation  on  the 
circuit  function  because  of  the  increasing  influence  of  tube  noise. 


APPENDIX  E 
MULTIPLE  SCATTERING  AND  FINITE  SCATTERER  CORRECTIONS 

The  slab  scatterers  used  in  the  experiment  were  between  l/5 
and  l/2  neutron  mean  free  paths  in  thickness.   Thus  a  significant 
number  of  multiply  scattered  neutrons  were  detected.  Cross  section 
data  must  be  corrected  for  this  effect.  A  second  correction  was  re- 
quired because  the  detection  crystals  were  close  enough  to  the  slab 
scatterer  to  detect  once- scattered  neutrons  over  an  angular  range  of 
about  ±8  degrees  in  polar  angle  with  respect  to  the  line  of  crystal- - 
scatterer  centers.   One  can  see  this  by  examining  Fig.  1  which  shows 
the  scatterer  and  detectors  in  typical  positions .   Because  the  detectors 
were  sampling  the  differential  cross  section  over  this  angular  spread, 
peaks  in  the  experimental  angular  distribution  were  reduced  and  valleys 
filled.  We  will  refer  to  the  correction  for  this  effect  as  a  "finite 
scatterer  correction." 

We  were  not  prepared  to  make  the  computer  analysis  necessary 
for  obtaining  accurate  corrections  because  for  cusp  analysis  a  first 
order  estimate  was  considered  satisfactory. 

The  measurements  of  differential  cross  sections  versus  energy 
were  made  for  only  one  or  two  angles  at  a  time.  Hence,  the  experimental 
points  were  correlated  in  energy  but  not  in  angle .  The  absolute  error 
in  the  (uncorrected)  measured  cross  section  was  approximately  ±15  per- 
cent (see  Error  Analysis)  and  was  independently  established  for  each 
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angle.  Therefore,  when  one  cross  plots  to  present  the  angular  distri- 
bution for  different  energies,  the  points  determining  the  shape  of  the 
angular  distribution  are  separately  uncertain  by  ±15  percent. 

The  shape  of  the  differential  scattering  cross  section  is,  of 
course,  essential  to  a  determination  of  the  multiple  scattering  correc- 
tion. We  used  first  order  multiple  scattering  and  finite  scatterer 

corrections  in  this  experiment  in  two  ways.  We  corrected  the  optical 

17 
model  prediction  of  Moldauer   for  the  angular  distribution  of  cerium 

and  compared  the  results  with  experimental  data  (see  Fig.  8).  We  also 

corrected  the  calculated  "energy-averaged"  cusp  effects  before  applying 

them  to  the  experimental  data  for  Ce,  Ba,  and  Zr. 

In  the  Fe  experiments,  the  structure  in  the  data  did  not  justify 
any  detailed  estimate  of  the  corrections,  short  of  a  full  computer  study 
on  the  angular  and  energy  dependence  of  multiple  scattering.  On  the 
basis  of  characteristic  scatterer  dimensions  we  estimated  an  upper 
limit  of  20  percent  for  the  corrections  in  the  case  of  iron. 

In  the  lithium  experiments,  the  scatterer  was  0.17  neutron 
mean  free  paths  in  characteristic  dimension.  This  dimension  is  small 
in  comparison  with  other  scatterer  thicknesses  typically  used  in  this 
experiment.  The  angular  distribution  is  peaked  backwards  and  the  lab 
neutron  energy  loss  for  back  scattering  is  from  30  to  Uo  percent  of 
incident  energy  (below  the  detection  bias  of  the  experiment).  All  of 
these  factors,  particularly  the  last,  contribute  to  an  estimate  that 
the  correction  for  multiple  scattering  in  this  case  is  small.  The 
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general  agreement  in  magnitude  between  the  measured  lithium  angular 
distribution  and  other  work  is  shown  in  Fig.  9  and  supports  this  con- 
clusion. 

It  remains  to  discuss  in  detail  the  first  order  multiple  scat- 
tering and  finite  scatterer  corrections  for  Ce,  Ba,  and  Zr.   We  used 

36 
a  procedure  and  notation  presented  by  Nauta.    Three  assumptions 

underlie  the  calculation:   first,  the  ratio  SK+1/S   is  independent 

of  K,   the  order  of  scattering,  where  S..  is  the  probability  of  scat- 
is. 

tering  elastically  K  times  before  leaving  the  scatterer;  second,  the 
angular  distribution  of  multiply  scattered  neutrons  is  the  same  as  in 
an  infinite  medium  of  the  same  material;  and  third,  there  is  no  loss 
of  lab  energy  upon  scattering.   The  scatterer  is  characterized  by  a 
parameter  6  which  is  the  fraction  of  neutrons  which,  having  been  once 
scattered  elastically,  have  undergone  at  least  one  more  interaction. 

Neglecting  third  and  higher  order  scattering,  the  necessary 

1   36 
formulas   are : 

da     =  f    (e)  da    ■  fl  +  (sJ^CeVs*1^))]'1  da   (0)   (El) 
corr    m.s.      exp         2  v  "    1  v  '  J  exp 

where 

00 

S2el(0)/Siel(0)  =  (2*6/oT  do  )  E  [2/(21  +  l)](a.)2  P.(cos  6) 

i=0 


1  -L 

a.  =  ^i  +  1)  /  dael  P.(u)du 


and  where  da      is  the  true  differential  cross  section;   da     is 
corr  '    exp 

the  measured  differential  cross  section;   a„  is  the  total  cross  section; 


13* 


and  da  .  =  a  ,  (Snel(0)  )/S-,el,  with  S-,el  =/  S,el(e)dft.   The  last 
el    el   1      '  1  1       1 

definition  emphasizes  the  fact  that  these  equations  are  circularly- 
related  and  that  approximations  of  some  kind  are  necessary  to  obtain 

a  solution.   Our  procedure  was  to  substitute  an  optical  model  angular 

17 
distribution   for  da  ;   thus  we  calculated  a  correction  to  the  opti- 
cal model  prediction  for  comparison  with  experiment.   Satisfactory 
agreement  between  the  corrected  optical  model  distributions  and  experi- 
ment was  considered  justification  for  applying  these  corrections  to 
the  comparison  of  predicted  cusp  effects  with  experiment.   This  general 
procedure  is  analogous  to  that  usually  adopted  in  solving  Eq.  (El); 
one  assumes  a  function  da  ,   and  by  iteration  of  Eq.  (El)  forces  it 

to  converge  to  da    .  We  stopped  this  process  at  a  first  iteration. 

corr 

In  essence,  the  method  separates  that  part  of  the  multiple  scat- 
tering which  depends  upon  the  angular  distribution  characteristic  of 
the  scattering  material  from  a  part  which  depends  upon  the  size,  shape, 
and  aspect  (with  respect  to  the  beam  axis)  of  the  scatterer.  The 
latter  dependence  is  contained  in  the  parameter  6. 

Lacking  a  better  estimate  of  5,  we  assumed  that  the  quantity 

cp  /cp    in  the  detailed  work  of  A.  Langsdorf,  Jr.,  R.  0.  Lane,  and 

37 

J.  E.  Monahan   is  1-5,  where  cp,   is  the  emergent  singly  scattered 

flux  per  unit  solid  angle  and  cp    is  the  unit  idealized  scattered 
flux  per  unit  solid  angle.   The  parameter  5  is  dependent  upon  the 
angle  of  scattering  at  which  the  cross  section  is  measured  because 
slab  scatterers  are  extremely  asymmetric.   In  the  work  of  Nauta,   who 
used  ring  scatterers,   6  could  be  assumed  independent  cf  angle. 
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Accordingly,  we  used  the  values  of  cp  /cp    given  in  Fig.  9  of 
Ref .  37  to  obtain  a  6  for  each  scatterer  and  for  each  angle  of  inter- 
est.  In  most  cases  extrapolation  of  the  curves  was  necessary.   These 
curves  are  entered  with  a  particular  value  of  h,   the  scatterer  thick- 
ness, in  units  of  neutron  mean  free  path. 

Our  scatterer  thicknesses,  calculated  at  the  neutron  energies 
of  interest,  were:   for  Zr,  H  =0.32;  for  Ce,  h  =  0.U6;  and  for  Ba, 

h  =  0.55- 

The  moments  a.   of  the  optical  model  angular  distribution  were 
included  to  a,  .  Multiple  scattering  corrections  as  a  function  of 
angle  were  then  computed  using  Eq.  (El)  and  the  appropriate  values  of 
6.   The  degree  to  which  corrected  optical  model  predictions  agreed 
with  experimental  data  is  indicated  in  Fig,  8.  We  estimated  the  over- 
all uncertainty  in  the  absolute  multiple  scattering  corrections  result- 
ing from  the  assumptions  and  approximations  in  the  foregoing  treatment 
to  be  ±30  percent.  Expressed  as  the  fraction  f    ,   which  is  con- 
venient for  modifying  data,  the  uncertainty  in  multiple  scattering 
corrections  adds  to  the  uncertainty  obtained  in  the  experimental  error 
analysis.  The  uncertainty  in  f      varies  with  scatterer  and  detec- 
tion angle,  as  indicated  in  Table  III. 

We  now  consider  the  calculation  of  finite  scatterer  corrections 
At  each  mean  polar  angle,  9   ,   once  scattered  neutrons  are  detected 
over  an  angle  spread  of  9     ±  A  9.      The  magnitude  of  2  A  9     was  de- 
termined graphically  from  the  geometry  of  the  experiment,  see  Fig.  1, 
and  was  found  to  vary  from  10  to  35  degrees.   The  angular  spread  is  a 
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function  of  9       because  the  scatterer-to-detector  distance  varied 
o 

depending  upon  the  experimental  situation  (counting  rate,  background  to 

effect  ratio)  and  also  because  the  projected  area  of  the  scatterer, 

normal  to  the  scatterer-detector  line  of  centers,  varies  with  9    . 

o 

As  in  the  correction  for  multiple  scattering,  the  procedure 
was  to  correct  an  optical  model  prediction  for  experimental  effects. 
The  finite  scatterer  correction  is  defined  by 

da     =  f _   (0)  da  (E2) 

corr    f.s.      exp 


f.p   (0)  "was  computed  using 


N 


f_   (e)  *  da  ./(l/N  S  da  +(sin  9   /sin  9   )  (E3) 

f.s.       opt'   '    .   opt      n      o 

n=l 

where  da  ^{9   )  is  the  optical  model  prediction  at  mean  angle  9   , 
opt  o  o 

and  the  denominator  represents  an  average  of  the  optical  model  distri- 
bution over  N  equal  divisions  of  the  spread  2  A  9     weighted  by  the 
solid  angle  subtended  by  the  nth  angle  division.   N  was  selected  to 
give  a  representative  average  and  varied  from  5  "to  8  depending  upon 
the  size  of  2  A  9. 

The  results  of  this  calculation  are  given  in  Table  III.  The 
magnitude  of  typical  finite  scatterer  corrections  is  indicated  on  the 
angular  distribution  shown  in  Fig.  8,  where  it  can  be  seen  that  the 
correction  is  small  except  in  the  valley  near  9      =  90° • 
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TABLE  III.  Multiple  Scattering  and  Finite  Scatterer  Corrections 


0 
o 

f 

m.s . 

ff.s. 

Uncertainty 

in  f    X  f 

m.s.   f .b. 

(percent) 

Cerium 

30° 

0.95 

1.0U 

3 

50° 

O.87 

0.98 

5 

70° 

0.80 

0.92 

8 

90° 

0.75 

0.90 

10 

110° 

0.78 

0.99 

9 

130° 

0.80 

1.00 

7 

150° 

O.85 

1.05 

6 

Barium 

90° 

O.65 

0.90 

18 

1250 

0.70 

1.00 

Ik 

Zirconium 

6o° 

0.90 

-1.0 

k 

130° 

O.85 

-1.0 

6 

Iron 

All 

0.8-1.0 

-1.0 

-10 

Lithium 

All 

-1.0 

1.0 

-  5 

The  factors  in  this  table  may  he  applied  to  the' j experimental 
cross  section  data  presented  in  this  report,  by  using 


da     =  f     f  _    da 
corr    m.s.   f.s.   exp 
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